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Abstract. We study the solution to the Robin boundary problem for 
the Laplacian in a Euclidean domain. We present some families of fractal 
domains where the infimum is greater than 0, and some other families 
of domains were it is equal to 0. We also give a new result on "trap 
domains" defined in [BCM], i.e., domains where reflecting Brownian 
motion takes a long time to reach the center of the domain. 

1. Introduction. 

The Robin problem (also known as the "third" boundary problem) for a Euclidean 
domain D C M. d is to find a function u such that 

A«(j) = 0, xeD, (1.1) 
dxi 

— = cu, x e dD, (1.2) 
on 

with one or more side conditions, where n is the unit inward normal vector field on dD, 
du/dn is the normal derivative of u in the distributional sense and c > is a constant. 
See Gustafson and Abe [GA] for the history of this problem. 

Our interest in the Robin problem stems from some recent applications in physics, 
electrochemistry, heterogeneous catalysis and physiology; see [FSF], [FS], [GFS], [Sa] and 
the references therein. Consider the mixed Dirichlet-Robin problem 

Au(x)=0, x £ D\B*, (1.3) 
On 

— — = cu, x e dD, (1.4) 
on 
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together with the side condition 

u{x) = 1, x G <9B*, (1.5) 

where 5* C D is a fixed closed ball with non-zero radius. The solution to (1.3)-(1.5) 
represents the steady state of a system in which some particles move randomly in D \ B* 
and cross a semi-permeable membrane dD. The other part of the boundary, dB*, is a 
source of particles and can be controlled so that we can assume a condition of type (1.5). 
The constant c in (1.4) is a physical characteristic of the membrane dD. One could consider 
a model with c dependent on x G dD but we will not do that in the present article. The 
constant c will play no role in our theorems so we will take c = 1 in the rest of the article. 

In some applied situations, it is desirable to have as much flux through the boundary 
as possible. The points of a man-made or natural membrane dD where there is no flux 
can be considered an inefficient use of material. Hence, it is interesting to know when the 
flux is non-negligible through all points of the membrane. In other words, we would like to 
know whether mf xe g D du/dn(x) > 0. In view of the relation (1.4) between the flux du/dn 
and the density u of particles and the maximum principle for the harmonic function u, 
this condition is equivalent to inf xeD \ Bf u(x) > 0. (By Lemma 2.4 below, we know u is 
non- negative.) 

Definition 1.1. We say that the whole surface of D is active if 

inf u(x) > 0. (1.6) 

xeD\B* 

If it is not the case that the whole surface is active, we say part of the surface is nearly 
inactive. 

In this paper we investigate the following problem. 

Problem 1.2. Give necessary and sufficient conditions of a geometric nature for the whole 
surface of D to be active. 

It is not difficult to show that the whole surface of a bounded Lipschitz domain is 
always active (see Remark 2.5(h) below). We have posed Problem 1.2 in terms of u rather 
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than du/dn because we are interested in non-Lipschitz domains D; so there are some 
boundary points where n is not well-defined while the solution u is always well-defined, 
and, in fact, is smooth in D \ B*. We do not have a complete solution to Problem 1.2, but 
we give a fairly explicit answer for some natural families of domains with fractal boundary. 

We will approach Problem 1.2 using probabilistic methods. This agrees well with the 
motivating physical models. Suppose that X is reflecting Brownian motion in D, L is 
its local time on dD, and T Br is the hitting time of by X. When D is a bounded 
C 3 -smooth domain, it is known that (see [MS] and [Pa]) 



This formula indicates that the third boundary problem (1.4) is more difficult to study 
from the probabilistic point of view than the corresponding Dirichlet and Neumann prob- 
lems. This is because the Dirichlet problem corresponds to killed Brownian motion and 
killing on the boundary presents no technical problems. The Neumann boundary problem 
corresponds to reflecting Brownian motion. The construction of reflecting Brownian mo- 
tion in an arbitrary domain D is a major technical challenge. Although this feat has been 
accomplished long time ago by Fukushima [Fu] on an abstract compactification, called 
the Martin-Kumarochi compactification, of D, many questions about the construction of 
reflecting Brownian motion on the Euclidean closure of a domain remain open (see [BBC]). 
Formula (1.7) shows that the Robin boundary problem (1.3)-(1.5) requires the construction 
and understanding of the local time. This is harder than constructing reflecting Brown- 
ian motion itself, because it is known that reflecting Brownian motion does not have a 
semimartingale decomposition in some domains. For some results in this area, see, e.g., 
DeB lassie and Toby [DT]. For information on the eigenvalue problem for the Laplacian 
with Robin boundary conditions, see Smits [Sml], [Sm2]. 

The following are the main contributions of this paper. The list includes some technical 
results that may have independent interest, 
(i) The solution of Problem 1.2 for a class of domains with fractal boundaries (Theorems 




(1.7) 



3.2 and 4.3). 
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(ii) A characterization of a class of "trap domains" in dimensions 3 and higher, improving 

a result in [BCM] (Theorem 5.1). 
(hi) Clarification of the rigorous meaning of solution to the differential equation (1.3)-(1.5), 
its existence, uniqueness, and probabilistic representation for non-smooth domains 
((2.3) and Lemma 2.4). In particular, we show that the solution to (1.3)-(1.5) is 
non-negative. 

(iv) A semimartingale decomposition of reflecting Brownian motion in a class of fractal 
domains (Theorem 2.2). 

(v) A sharp estimate for the Green function with Neumann boundary conditions in long 
and thin domains (Lemma 4.4). 

(vi) A new version of the Neumann boundary Harnack principle, stronger than the one in 
[BH] (Lemma 2.8). 

(vii) The proof that reflecting Brownian motion starting from the cusp point is not a semi- 
martingale, for some cusps (Remark 4.14). This complements a result of Fukushima 
and Tomisaki [FT]. 

A simple example illustrating our main theorems is a cusp domain, defined for a fixed 
a > 1 by 

D = |x = (x!,X2, ■ ■ ■ ,Xd) ■ < xi < 1 and xf > {x\ + . . . + a^) 1/2 j • 

Applying the main results (Theorems 3.2 and 4.3) of this paper, we show in Example 
3.4 (for d = 2) and Example 4.13 (for d > 3) that the whole boundary of D is active if 
a E (1, 2), and part of dD is nearly inactive if a > 2. There are more examples given in 
Sections 3 and 4. 

The paper is organized as follows. Section 2 contains some technical preliminaries, 
many of which may have independent interest. Section 3 presents the solution to Problem 
1.2 for a class of 2-dimensional domains, using techniques developed in [BCM]. Section 
4 is devoted to Problem 1.2 in dimensions 3 and higher. Finally, Section 5 presents an 
application of the techniques developed in Section 4 to "trap domains" in dimensions 3 
and higher. 
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2. Reflecting Brownian motion in domains with fractal boundaries and the 
Neumann boundary Harnack principle. 

This section is devoted to two important technical aspects of this paper. First, we 
will show that reflecting Brownian motion has a semimartingale decomposition for a class 
of fractal domains that contains some natural examples. The second technical result is a 
boundary Harnack principle for harmonic functions satisfying Neumann boundary condi- 
tions. 

We will let | • | stand for the Euclidean norm in R. d (for any dimension d > 1), for the 
volume ((i-dimensional Lebesgue measure) of a set A C M d , and for the (d — l)-dimensional 
surface area of the boundary d A of a set A C M d . The meaning will be obvious from the 
context so this notation should not lead to any confusion. For an open set D of R d , C C (D) 
and C^°(D) denote the space of continuous functions with compact support in D and the 
space of smooth functions with compact support in D, respectively. 

A ball with center x and radius r will be denoted B{x, r). The notation will refer to 
an open ball, unless noted otherwise. 

The harmonic measure of a set A C dD in the domain D, relative to z, will be denoted 
u(z, A, D). 

The distribution of Brownian motion in D \ B* starting from x G D \ reflected on 
dD, and killed on dB* will be denoted P x . The corresponding expectation will be denoted 
E x . The hitting time of a set A will be denoted T A , i.e., T A = inf{£ > : X t G A}. We 
will sometimes write T* or Tjf to show the dependence of the hitting time on the process. 

We will use elements of excursion theory and Doob's /i-processes. See [D] for the 
discussion of /i-transforms in the case of (non-reflecting) Brownian motion, and [Sh] for 
conditioning of general Markov processes. Elements of excursion theory can be found in 
[Mv], [Bl], [Bu] and [Sh]. 

A real-valued function / defined on A C R d is called Lipschitz with constant A < oo 
if \f(x) — f(y)\ < X\x — y\ for all x,y G A. A domain D is called Lipschitz if there exist 
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r > and A < oo such that for every x G 3D, the set dD n B(x,r) is the graph of a 
Lipschitz function with constant A in some orthonormal coordinate system. We call (A, r) 
the Lipschitz characteristics of D. 

Definition 2.1. We will say that a domain D belongs to class T> if there exists an increasing 
sequence of domains D n C D with the following properties. 

(i) Each D n is a Lipschitz domain with characteristics (A, r n ) (all the A's are the same 
but the r n 's may differ) and U^Li D n = D. 

(ii) For every n > 1, the set <9.D n D dD is a subset of the relative interior of dD n+1 fl dD. 

(iii) sup n>1 \dD n \ < oo and lim,^^ \dD n \ dD\ = 0. 

The set d^D = \J n dD n fl dD will be called the Lipschitz part of dD. 

Every bounded Lipschitz domain is in V. See Examples 3.4, 3.6, 4.13 and 4.14 below 
for domains D ET> which are not Lipschitz. 

Constructing a reflecting Brownian motion on a non-smooth domain D is a delicate 
problem. Let 

W 1,2 {D) = {/ G L 2 (D,dx) : V/ G L 2 (D,dx)} 

be the Sobolev space on D of order (1,2). Fukushima [Fu] used the Martin-Kuramochi 
compactification D* of D to construct a continuous diffusion process X* on D* with 
transition semigroup denoted Pt, such that 

{/GL 2 (D,^):sup^ / f(x)(f(x)-P t f(x))dx<oo} = W 1 ' 2 (D) 
t>o t Jd 

and for / G W 1 - 2 ^), 

£(/,/) = limi f f{x){f{x)-P t f{x))dx=\ [ \Vf(x)\ 2 dx. 

The pair (5, W^ 2 (D)) is called the Dirichlet space ofX* in L 2 (D*,m), where m is Lebesgue 
measure on D extended to D* by setting m(D* \ D) = 0. See [FOT] for definitions and 
properties of Dirichlet spaces, including the notions of quasi-everywhere, quasi-continuous, 
etc. The process X* could be called reflecting Brownian motion in D but it lives on an 
abstract space D* that contains D as a dense open set. Chen [CI] proposed referring to the 
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quasi-continuous projection X of X* from D* into the Euclidean closure D as reflecting 
Brownian motion in D. The projection process X is a continuous process on D, but in 
general X is not a strong Markov process on D (for example this is the case when D is the 
unit disk with a slit removed). However when D is a Lipschitz domain, it is shown that X 
is the usual reflecting Brownian motion in D as constructed in [BH]. It was proved in [CI] 
that, roughly speaking, if dD has "finite surface measure," then X is a semimartingale 
and has a Skorokhod decomposition. This result was further sharpened in [CFW]. See 
the introductions of [CI] and [CFW] for the history of constructing reflecting Brownian 
motion on non-smooth domains. 

Theorem 2.2. If D E V, then reflecting Brownian motion X in D starting from x E 
DUOlD has a semimartingale decomposition X t = x+Wt+N t , where Wt is a d-dimensional 
Brownian motion, 

N t = I n(X s )dL s , 
Jo 

and L, the local time, is a non-decreasing continuous process that does not increase when 
X is not in OlD, i.e., J Q l^Q L D) c (X t )dL t = 0. The Revuz measure of L for the process 
X* is surface measure on 8lD. 

Note that the local time L in our theorem satisfies the condition J °° l^g L ^ c (Xt)dLt = 
0, which is stronger than the usual condition J °° lr){X t )dL t = 0. 

Proof. Let {D n , n > 1} be the increasing sequence of Lipschitz domains in the definition 
of D E V. Let D* be the Martin-Kuramochi compactification of D used in [Fu]. To be 
precise, for every a > 0, let H a denote the space of all h in D such that (a — \&)h = in 
D and having 

£ a {h,h) = - j \Ah(x)\ 2 dx + f u(x) 2 dx < oo. 
2 J d J D 

For y E D, let x i— > H a (x,y) be the unique a-harmonic function in H a such that 

£ a {H a { ■ , y), v(-)) = v(y) for every v E H a . 

Let G° a (x,y) be the ct-resolvent density function for Brownian motion in D killed upon 
exiting D. Define 

G a (x, y) = G° a (x, y) + H a (x, y). 
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It is shown in [Fu] that x h- > G a (x, y) is continuous on D \ {y} and G a (x, y) = G a (y,x). 
Define a metric 5 on D by 

S(x,y)= [ (\G 1 (x,z)-G 1 (y,z)\Al)dz 

and let D* be the completion of D under the metric S. Fukushima [Fu] showed that there 
is a conservative continuous Hunt process X* on D* \ N associated with the Dirichlet 
space (£, W 1,2 (D)) on L 2 (D*,m), where N is a set that has zero capacity with respect to 
(£, W 1,2 (D)) and m is Lebesgue measure on D extended to D* by defining m(D* \ D) = 0. 
Since each coordinate function Xi G W 1,2 (D), then each coordinate function admits a quasi- 
continuous version on D* , which will be denoted as fi. Note that (/i,- • •, fd) is defined 
quasi-everywhere on D* and is a quasi-continuous map from D* into D. Define 

X=(f 1 (X*),...,f d (X*)). 

Then X is a conservative continuous process on D, which is called reflecting Brownian 
motion on D in [CI]. It coincides with the usual reflecting Brownian motion when D is a 
bounded Lipschitz domain. 

Let X n be reflecting Brownian motion on D n . It is known from [BH] that X n has 
a Holder continuous transition density function p n (t,x,y) on (0, oo) x D n x D n . Its a- 
resolvent density function will be denoted as G n (x,y). Define 

r n = mf{t > : X™ G <9D n \ <9L>} 

and 

G^°(x,y ) = G a (x,y ) - K [e- aT »G a (X? n ,yo)] . 

It is easy to verify that is the ct-resolvent for reflecting Brownian motion in D n killed 
upon hitting dD n \ dD. Thus we have 

G a (x, yo) = GS(x, ?/) + E£ [e"^" (G tt W B , W>) - G^X^yo))} . 

By [BH], x i— > J/o) is continuous on D n \ {y } and 

x ^ E£ [ e — (G«(X» , y ) - G2(X» , y ))] 
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is continuous on dD n fl dD since it is harmonic in D n \ {y } with zero Neumann boundary 
conditions on dD n n <9-D and zero Dirichlet boundary conditions on <9D n \ dD. Hence 
we conclude that x h- > G(x,yo) extends continuously to <9-D n fl <9.D under the Euclidean 
topology for every n > 1 and hence to <9l-D = {J^ =1 (dD n fl <9.D). This implies that 

d L D C (L>* \ D) n <9D. 

Note that on <9l-D, (/i, • • • , /d) is the identity map and so = X t when X t * £ 8lD. 

Let ct/j denote surface measure on dD k and let n k (x) be the unit inward normal vector 
field on dD k which is defined almost everywhere with respect to a k . By the definition of 
D £ V, 

k i-> a k (dD k ndD) 

is an increasing function and 

lim a k (8D k ) = lim a k (8D k n dD) = a(8 L D), 

since dD k ndD C dD k+ iHdD and lim^oo a k (dD k \dD) = 0. Here a is surface measure on 
<9l_D. Since sup fc>1 a k (dD k ) < 00, there exist a subsequence {fcj, J > 1} and finite signed 
measures (z/ 1? • • • , z/ d ) on D* such that n fcj .cr fcj converges weakly on D* to (z/ 1? • • • , u d ); that 
is, 

lim / (gi(x), ■ ■ -,gd(x)) ■ n k j (x)a k Adx) = V] / ^(x) z^(cfo), (2.1) 

for all bounded continuous functions {<7i, • • • , on D* . For every 1 < i < d and k > 1, 
H(Z>*\(d£> fc ndZ>)) < lim a kj (D*\(dD k ndD)) 

= lim (T fc .(aD fc , \(aD fc na^)) 

= <r(d L D\dD k ). 

Thus for 1 < z < d, 

\vi\{D* \ d L D) = lim H(D* \ (d£> fc n dD)) < lim ct(<9 Lj D \ dD k ) = 0. (2.2) 

A; — >oo k — >oo 

On the other hand, by the definition of D £ £>, n k a k converges weakly on D to ncr, where 
n is the unit inward normal vector field of D on 8lD in the following sense: 
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k^ooJ D j D 



lim /_(^i(a;), • • • , g d (x)) ■ n k (x) a k (dx) = /_(#i(ir), • • • , </d(x)) • n(x) a(dx) 
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for all bounded continuous functions {gi, • • • , go] on D that vanish on dD n \ 3D for some 
n > 1. 

Since d L D cD*n5, we conclude from (2.1) and (2.2) that 

(fi, • • • ,^d) = ncr on D*. 

By Theorem 4.4 of [CI], a is a smooth measure of X* and thus it determines a positive 
continuous additive function L of X*. Moreover, 

X t =X + Wt+ [ n(X s )dL s for t > 0, 
Jo 

where W is a <i-dimensional Brownian motion. The above Skorokhod decomposition holds 
for quasi-every starting point Xq in D* with Xq = f(Xo). However, since the a-resolvent 
density function x i— > G(x,yo) is continuous on d^D U (D \ {yo}), reflecting Brownian 
motion X* can be defined to start from every point x G D fl 8lD (cf. [FOT]). Hence the 
above Skorokhod decomposition holds for every starting point Xq E D Li 8lD. Clearly, 
since a is carried on <9l-D, 

/ 1 {x s ^d L D}dL s = l{x;$d L D}dL a = Q. 
Jo Jo 

This proves the theorem. □ 



Remark 2.3. Let td\j8 l d be the first exit time from D U &lD by reflecting Brownian 
motion on D. Starting from x G D U <9l-D, {X t ,t < t^hj^d} is a strong Markov process 
onDU <9l-D, since it coincides with {X?, t < tduo l d}- Here 

r D ud L D = mf{t >0:X t (£DU 8 L D} = mf{t >0:X?(£DU 8 L D}. 

However even under the conditions of Theorem 2.2, reflecting Brownian motion on D 
may not be a strong Markov process. For example, let D be the union of {(x,y) G M 2 : 
\y\ > \x\ and \y\ < 1} and {(x, y) G M 2 : 1 < x 2 + y 2 < 4}. Then clearly reflecting 
Brownian motion X on D can not have the strong Markov property since when X t is at 
the origin 0, one can not tell how it will be reflected unless one knows where it came from. 
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Of course, in this example, for starting points in D other than the origin 0, reflecting 
Brownian motion will not visit 0. But one can modify this example so that the set of 
such non-Markovian points has positive capacity so it will be visited by the reflecting 
Brownian motion. Here is such an example. Let K be the standard Cantor set in [0, 1]. 
Let C = {(x, y) G M 2 : |x| < \y\ < 1}. Define D C M 2 by 

D = {(x,y):Kx 2 + y 2 < 9} U |J ((x,0) + C)). 

xeK 

Clearly D satisfies the assumptions of Theorem 2.2. Note K = K x {0} is the set of non- 
Markovian points and K has positive capacity (see [C2]) and so will be visited by reflecting 
Brownian motion in D. □ 



Now we make precise the meaning of solution to the partial differential equation with 
Robin and Dirichlet boundary conditions (1.3)-(1.5) for D G V. Define 

W 1,2 (D; B m ) = {u G W X ' 2 {D) : u = q.e. on the closed ball B*}. 

We say u is a (weak) solution of (1.3)-(1.5) if the following two conditions are satisfied. 

(i) u and its distributional derivative Vu are in L 2 (D \ B*) and for any bounded g G 
W^iD-B*), 

/ V g(x) -Vu{x) dx = —c / g(x)u(x)a(dx). (2.3) 

JD\B, Jd L D 

(ii) u is continuous in a neighborhood of dB* and u = 1 on dB*. 

Note that any / G W 1,2 (D; B*) admits a quasi-continuous version on D* \ B*. 
Throughout this paper, we will always represent such / by its quasi-continuous version, 
which will still be denoted as /. In particular, / is well defined q.e. on D* \ D. Since 
8lD C D* \ D and a is a smooth measure of X* according to Theorem 2.2, / is well 
defined cr-a.e. on 8lD for every / G W 1,2 (D; £?*). Hence the right hand side of (2.3) is 
well defined. 
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Lemma 2.4. The partial differential equation with Robin and Dirichlet boundary condi- 
tions (1.3)-(1.5), where c > is a constant, has a unique solution u(x) given by u(x) = 
E x [exp(— %Lt b * )] • I 11 particular, u is non-negative. 

Proof. We first establish existence. Note that u(x) = E x [exp( — §Lt b „ )] is well defined for 
every x G D U 8lD and for q.e. x G D* . By the Markov property of X* , for iGDU d^D 
(as well as for q.e. other x G D*), 



v(x) = 1 - ~E X exp(-|L Ts J 



— — E T 
2 



Let X*'° be reflecting Brownian motion X* killed upon hitting B* and let the transition 
semigroup be denoted by {P t °,t > 0}. It is known (cf. [FOT]) that the Dirichlet form of 
X*<° is (S, W^ 2 (D; £*)) on L 2 (D* \ B*, m). For q.e. x G D* \ £*, 



- P t °v(x) = -E a 



Hence 



lim - / v(x)(v(x) — P?v(x)) dx = — lim / w(x)E x 



u(X s *'°)c2L 



= - / w(a;)w(x)(j((ix) < oo. 
2 ^-D 



Thus v G VF 1,2 (.D; £?*), and a similar calculation to the above yields that for any bounded 



z •/ D\B* 



Vg(x) ■ Vv(x) dx = £(g, v) 



lim - 



t Jd\b« 
°- \ g{x)u{x)o{dx) 

1 Jd L D 



g(x)(v(x) — P®v(x)) dx 



'Si J 

This shows that v is harmonic in D \ B* and = — eu. In particular, i> is continuous 
in D\B*. Since every point of d-B* is regular, we see that v vanishes continuously on 
dB*. Translating these properties to the function u = 1 — v shows that u is a solution to 
(1.3)-(1.5). 
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Now we show the uniqueness. Suppose that u\ and u 2 are two solutions for (1.3)-(1.5). 
Define w = u\ — u<z- Then w e W 1 ' 2 {D; B*) and it follows from (2.3) that 

/ Vg{x) ■ Vw{x) dx = — c / g(x)w(x)a(dx). 
Jd\b* Jd L D 

Letting g = ((— n) V w) A n in the above and then letting n — > oo, we have 



' D\B 

Since c > 0, we must have 



/ | Vw(x) | 2 (ix = — c / 

Jd\B* Jd L D 



w(x) \ 2 cr(dx). 



[ \Vw{x)\ 2 dx= f \w(x)\ 2 a(dx) = 0. (2.4) 

•/ D\B* Jd L D 

Since D \ B* is connected, to has to be constant in D \ , while the second equality in 
(2.4) implies that w = cr-a.e. on Therefore w = in D \ B* and hence u\ = u 2 . 

This establishes the uniqueness and completes the proof of this Lemma. □ 



Remark 2.5. (i) A simple modification of the above argument establishes the existence 
and uniqueness for solutions to the Robin problem (1.3)-(1.5) with c being a bounded 
non-negative function. The solution in this case can be represented as 



u(x) = E a 



exp 



c(X s )dL t 



for x G D \ B*. 



(ii) Suppose that D is a bounded Lipschitz domain in R d with d > 3 and u is the 
solution to (1.3)-(1.5). By Jensen's inequality, we have 



u(x) > exp ^— -E x Lt Ss( j for x E D \ B* 



Let Gd\b* be the Green function of the reflecting Brownian motion killed upon hitting 
5*. It is known from [BHl that 



G d \bX x ,v) < 



\x-y 



d-2 



for x, y e D \ B*. 
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It follows then 



sup E x L TBtt = sup / G D \ B ,(x,y)a(dy)< sup / Cl ld _ 2 a(dy) < oo. 

xeD\B* xED\B*JdD xeD\B„ J dD \ x ~ V\ 

Hence ini xeD \ Bt u(x) > 0. In other words, the whole surface of a bounded Lipschitz 
domain in R d with d > 3 is always active. □ 



Let D C l d be a Lipschitz domain and let O be a connected open set in R d . The 
following definition of "Neumann boundary conditions" for a harmonic function is standard 
in analysis and PDE (cf. [K]). 

Definition 2.6. A function h defined on D n O is said to be harmonic in D n O with 
zero Neumann boundary conditions on dD n O if ft G W 1,2 (Oi H D) for every relatively 
compact open subset Oi of O and 

/ V/i(x) • Vip(x)dx = (2.5) 

for every ip G C%°(Oi) and consequently for every continuous ip G VF 1,2 (Oi) that vanishes 
on dO\. 

The following lemma says that functions expressed in terms of the hitting distribution 
of reflecting Brownian motion in D are harmonic functions with zero Neumann boundary 
conditions in the sense of Definition 2.6. 

Lemma 2.7. Let X be reflecting Brownian motion in the Lipschitz domain D, and O a 
connected open subset ofR d . Define r = mf{t > : X t D n O}. Then for any bounded 
measurable function ip on dO fl D, 

h(x) = E x [iJj(X T )], xeDHO, 

is a harmonic function in D fl O with zero Neumann boundary conditions at 3D fl O. 

Proof. Without loss of generality, we may assume that ip > 0. Define 

df / X t if t < t 

1 \d iit>r, 
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which is reflecting Brownian motion in D killed upon leaving O. It is well-known that X° 
is a symmetric Markov process on D n O with Dirichlet form (£, VF 1,2 (-D; O c )), where 

W 1 ' 2 ^; O c ) = {u G W 1 ' 2 ^) : u = q.e. on O c }. 

The transition semigroup for X° will be denoted by {P® ,t > 0}. 

Let Oi be a relatively compact open subset of O and let / > be C\ with supp[/] C O 
and / = 1 on 0\. Define u(x) = f(x)h(x). Then for x G D n O, 

u(a;) - P t °w(x) = E x [(/(X ) - f(X t ))h(X t ); t < r] + E x [f(X )h(X T ); t > r] . 

Note that by time-reversal, 

/ u(x)E x [(f(X ) - f(X t ))h(X t ); t<r]dx 
JDno 

= [ E x [f(X )h(X )(f(X ) - f(X t ))h(X t ); t<r]dx 
J DnO 

= [ E x [f(X t )h(X t )(f(X t ) - f(X ))h(X ); t < r] dx. 

J DnO 



Hence 



'DnO 

1 

2 



DnO 

2 



Thus 



/ u(x)E x [(f(X ) - f(X t ))h(X t ); t<r]dx 
JDno 

E x [(f(X ) - f(X t )) 2 h(X )h(X t ); t<r]dx 

<Mi / E x [(f(X t ) - f(X )) 2 ; t < r] dx. 
1 JDno 

limsup- / u(x)(u(x) — P^u(x))dx 
t— o t J Dn0 

< i imsup (HI f Ex [(f(X t ) - f(X )) 2 ; t<r]dx 
t— o v It J Dn o 

f(x) 2 P x (t > r)dx) 



DnO 



< \\h\\lo 



[ \Vf{x)\ 2 dx < oo, 
J DnO 



'DnO 

by Lemma 4.5.2(i) and (4.5.7) of [FOT]. This implies, by Lemma 1.3.4 of [FOT], that 
u G W 1 ' 2 (D;O c ) and so h G W 1 ' 2 (Oi). 

u(x) = E x u(X° ) for x G O n D, 
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where r 0l = inf{* > : X® £ Oi C\ D}. Hence by Theorem 4.3.2 of [FOT], u is S- 
orthogonal to W 1 ' 2 (D; Of); that is, 

\f Vu(x) -V(t>{x)dx = Q for every e W ia {D;0{). 

This shows that u and therefore h is harmonic in D H 0\ with zero Neumann boundary 
conditions on dD fl 0\. Since Oi is an arbitrary relatively compact open subset of O, we 
conclude that u is harmonic in DnO with zero Neumann boundary conditions on dDDO. 
□ 

The following version of the Neumann boundary Harnack principle is similar to (but 
slightly more general) than Theorem 3.9 of [BH]. The result in [BH] is limited to smooth 
domains whose boundaries are locally graphs of Lipschitz functions (although the constant 
in that theorem depends only on the Lipschitz constant A) and to harmonic functions h as 
in our Lemma 2.7, with non-negative ip. 

Lemma 2.8 (Neumann boundary Harnack principle). Suppose that $ : IR d_1 — > R 
is a Lipschitz function with constant X < oo, i.e., \$>(x)— <&(y)| < \\x—y\ forallx,y G 
Assume that $(0) = and let D = {x = (x\, x 2 , • • • , G R d : £d > $>((xi, . . . , Xd_i))}. 
If r > 0, C\ > 1, and h : 5(0, Cir) fl D — > [0, oo) is harmonic with zero Neumann boundary 
conditions on £?(0, Cir) fl <9.D tiien 

h(x)>c 2 h(y) for all x,y G B(0,r) HD, (2.6) 

where C2 > depends only on A and c\. 

Proof. For (j/i, • • • , j/d) G lR d , denote y = (yi, • • • , yd-i)- Define a one-to-one map 4> : 
(j)(y,yd) = (y>y<i ~~ ^(y))- As ^ * s Lipschitz, the Jacobians of (p and its inverse </> _1 are 
bounded, with the bound depending only on the Lipschitz constant A. Under 0, |A is 
mapped into a uniformly elliptic divergence form operator L with coefficient matrix A(x) 
(see Remark 2.1.4 of [K]). Let U = 0(5(0,0^) n 75) and u{x) = h{<fr\x)) for x E U. 
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Using the change of variable formula, we conclude from (2.5) that for every continuous 
ip G W X ' 2 {U) that vanishes on dU n {y G R d : y d > 0}, 

A(x)Vu(x) -Vi)(x)dx = 0. (2.7) 

Let U~ be the "mirror" reflection of U with respect to the hyperplane {(y, yd) '■ yd = 0}, 
that is, U~ = {y = (y,y d ) : (y, -yd) e F °r 2/ = (y,Vd) e define A(y) = 

A((y, — y^)) and it(y) = u((y, —ya))- Then L = V(AV) is the uniformly elliptic diver- 
gence form operator defined on the domain U U U~ . It now follows from (2.7) and its 
corresponding version for U~ that 

I A(x)Vu(x) ■ Vip(x)dx = for every V G C™(U U U~). 
Juuu- 

Hence u is a non-negative L-harmonic function on UUU~ . The desired Harnack inequality 
for h now follows from the Harnack inequality for the L-harmonic function u. □ 



Remarks 2.9. (i) Some regularity conditions for a harmonic function with zero Neumann 
boundary conditions have to be assumed (such as those formulated in Definition 2.6) in 
order for the Neumann boundary Harnack principle to hold, even if D has a C°° boundary. 
The Neumann boundary Harnack principle does not need to hold for a harmonic function in 
DnB(io, r) which satisfies zero Neumann boundary conditions only almost everywhere on 
dD fl B(xo,r). For example, let D be a half-space in IR d , d > 3, with dD passing through 
the origin, and let h(x) = \x\ 2 ~ d . Then h satisfies the Neumann boundary conditions 
everywhere except at the origin. The Neumann boundary Harnack principle does not hold 
for this function h in D fl B(0, 1). 

(ii) We will apply Lemma 2.8 to two classes of functions. One of these families consists 
of harmonic functions defined in a probabilistic way, as in Lemma 2.7. That lemma shows 
that Lemma 2.8 is applicable to harmonic functions in this family. 

We will also apply Lemma 2.8 to the Green function x — > G(x, y), where y G D\ B*, 
and G{ • , y) is the density of the expected occupation measure for the reflecting Brownian 
motion in a Lipschitz domain D killed upon hitting starting from y. To see that 
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Lemma 2.8 can be applied, consider any y E D \ and let U be any relatively compact 
subdomain of D \ (B* U {y}). Then for x E U, G(x,y) = ~E x [G(X Tu , y)]. So by Lemma 
2.7, x — > G(x,y) is "locally" in W 1,2 (D) and is harmonic with zero Neumann boundary 
conditions on dD. 

3. Simply connected planar domains. 

This section will present some results based on ideas developed in [BCM], a paper 
on "trap" domains. We will present a new result on trap domains in Section 5. In this 
section, we will review only as much of the material from [BCM] as is relevant to Problem 
1.2. We will use complex analytic notation and concepts. Consult [Po] for the definitions 
of prime ends, harmonic measure, etc. 

We start with some definitions that apply to domains in any number of dimensions. 
Let X be normally reflecting Brownian motion on D C K d , d > 2, starting from x E D 
and killed upon hitting a closed ball £?*. As is mentioned in the previous section, X is 
obtained as the projection of reflecting Brownian motion X* on the Martin-Kuramochi 
compactification of D into D. The distributions of both X and X* will be denoted P x 
and the corresponding expectations will be denoted E x . Let G{x,y) be defined on (D \ 



where dy denotes (/-dimensional Lebesgue measure. Clearly G(x, y) is a symmetric function 
on (D \ B*) x (D \ B*). It follows from Lemma 3.2 of [CFW] that the function G(x,y) 
can be extended continuously to (D* \ B#) x (D\B*), where D* is the Martin-Kuramochi 
compactification of D as mentioned in the proof of Theorem 2.2 in the previous section. 
Note that (cf. Section 2.1 of [BCM]) if x, y E D \ B* and x ^ y then t i-> G(X*, y) is a 
continuous local martingale. It is easy to see that G(x, y) is the Green function for the 
domain D\B* with (zero) Neumann boundary conditions on dD (in the distributional 
sense) and (zero) Dirichlet boundary conditions on dB*. 

For the rest of this section, suppose that D is a simply connected open subset of the 
complex plane C, z* is the center of and zq is a prime end in D. Consider a collection 



B*) x (D \ B*) by 
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{ln}n>i of non-intersecting cross cuts of D that do not intersect and such that 7 n +i 
separates 7 n from z$ and the 7 n 's tend to Zq. Suppose further that a is a curve in D 
connecting z* to zq such that a fl 7 n is a single point z n , for each n. This system of curves 
divides D into subregions: let Q n denote the component of D \ 7 n which does not contain 
z*. Thus D n = fl n \ O n _|_i is the region between 7 n and 7 n +i- Write Oi \ a = fi + U Q~ , 
where each set Q + and fi~ is connected, and set -D+ = + fl D n and _D~ = fi~ fl D n . See 
Figure 3.1. 




Figure 3.1. Hyperbolic blocks. 



Recall that the harmonic measure of a set A C <9-D in the domain D, relative to z, is 
denoted u(z, A, D). 



Definition 3.1 We will say that the system of curves {7 n } U a divide D into hyperbolic 
blocks tending to the prime end z$ if for some c* > and all n > 1, the following conditions 
hold: 

(i) c* < v{z*, dfl+ n dD, D) < 1/2 and c* < 90" n d£>, D) < 1/2, 

(ii) for all n > 1 and for all z G <9-D^ U {,z n _i}, we have a; (2;, <9-D^ fl dD, D) > c*, 
(hi) for all n > 1 and for all z G U {z^-i}, we have u(z, dD~ fl dD, D) > c*. 

We will call a system of hyperbolic blocks regular if it satisfies in addition the following 
condition, 
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(iv) for every n > 1 there exists z G 8D+ n <9-D n such that u;(2;, <9.D n D <9-D, _D n ) > c* and 



For every simply connected domain and any prime end zq, there exists a family of 
regular hyperbolic blocks. Here is one way to construct {7 n }n>i and a. Suppose that <p 
is a conformal map of the upper half plane H onto D, such that f(0) = zq and ip(i) = z*. 
Then we can take j n = <p(B. n {\z\ = 2" n }), n > 1, and cr = {^(zy) : < y < 1}. The 
conformal invariance of harmonic measure makes it is easy to verify that the ensemble 
{^n} U a divides D into hyperbolic blocks tending to zq. Condition (iv) is satisfied by 
z = </?(i(3/4)2~ n ). Hyperbolic blocks are useful because they can be constructed geomet- 
rically, without knowledge of any properties of the mapping <p; see [BCM] for examples of 
hyperbolic blocks. 

In typical examples, verifying conditions (i)-(iv) is not harder than verifying just (i)- 
(iii). We did not include (iv) in the definition of hyperbolic blocks in order to keep the 
same nomenclature as that in [BCM]. 

Theorem 3.2. Let D e V be a simply connected planar domain, 
(i) If there exist constants c*,c' G (0, oo) such that for each prime end z G dD there is 
a system of curves {~y n } U a dividing D into regular hyperbolic blocks with parameter 



then the whole surface of D is active, 
(ii) Let r n denote the distance between 7 n and 7 n +i- If for some prime end zq G dD, 
there is a system of curves {^ n } U o dividing D into regular hyperbolic blocks with 



u(z,8D+ f]8D,D n ) > c*. 



c* and 




n 



(3.1) 



oo 




(3.2) 



n=l 



then part of the surface of D is nearly inactive. 



Example 3.4 and especially Example 3.6 show that the gap between parts (i) and (ii) 
of Theorem 3.2 is not large. 
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We need a lemma to prove Theorem 3.2. 



Lemma 3.3 Suppose that D E V and the {7^} divide D into regular hyperbolic blocks. 
For n > 1 and y E 7 n -i, let x 1— > be tie Poisson kernel with pole at y for reflecting 

Brownian motion in f2 n _i killed upon hitting 7 n _i. Let P% denote the distribution of 
Doob's hy-transform of reflecting Brownian motion on D* killed upon hitting 7 n -i 3 starting 
from x E Note that P|-a.s., the process will stay in the closure of O n _i in D* until 

its lifetime. Let r n denote the distance between 7 n _i and 7 n . There exist ci,pi > 0, 
depending only on D and c*, such that P^(Ly i > Cir n ) > pi for any x G 7 n . 

Proof. Let </? be a one-to-one conformal map of -D n _i onto the unit disc S 1 = G 
C : |z| < 1}, such that V?(dL>+-i n dL>) = 7i = {z = e l6 : 9 X < 9 < n - 0i} and 
(p{dD-_ x r\dD) = I 2 = {z = e~ ld :9 1 <9 < 7r-0i}, for some < 9 1 < tt/2. By condition 
(iv) in Definition 3.1 and conformal invariance, there exists a point xi E S such that 
uj(xi, ii, S) > c* and u;(xi, I 2 , S) > c*. This easily implies that there exists 9 2 = 6> 2 (c*) G 
(7r/4,7r/2) such that X < 2 - Let 3 = (tt/2 + 2 )/2, Ji = {2 = e* e : 3 < < ^ ~ 9 3 } 
and J 2 = {z = e~ ld : 3 < 9 < n — 3 }. For some c 2 = c 2 (c*) > and every z G S on 
the imaginary axis, cj(z, Ji U J 2 , S) > c 2 . Let J r and J be the right and left connected 
components of dS \ (I\ U 7 2 ). Let X be Brownian motion in S 1 with normal reflection on 
1 1 U I 2 killed upon hitting J r Li J e . It is easy to see that for some c 3 = 03(0*) > and 
every z G J\ U J 2 , if reflecting Brownian motion in £ starts from z, then it hits J £ before 
hitting J r with probability greater than c 3 but less than 1 — c 3 . 

Let 7 n _i/2 = ^P~ 1 ({z = a + bi E S : a = 0}), K\ = </? -1 (Ji) and K 2 = </? _1 (J 2 ). By 
conformal invariance, for every x E 7 n _i/ 2 , we have uj(z, K\ U K 2 ,D n _i) > c 2 , and for 
every point x E K\ U K 2 , the probability that reflecting Brownian motion in D starting 
from x hits 7 n _i before hitting 7 n is in the range (c 3 , 1 — c 3 ). 

Find c 4 G (0, 1/8) so small that a Brownian motion W starting from x will make a 
double loop in an annulus B(x,r) \ B(x,3r/4) for some r E (c±r n , r n /8), and then will 
make a crossing from B(x + zr/2,r/16) to the ball B(x + i2r, r/16) within the convex 
hull of the two balls, all before leaving B(x, r n /3), with probability greater than 1 — c 3 /2. 
A "double loop" in B(x,r) \ B{x, 3r/4) means that there exist t\ < t 2 such that W t E 
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B(x,r) \ 3r/4) for all t G (£1,^2), and a continuous version of t — > arg(W"t — x) 
increases by 4-7T over the interval [£1, £2]- Note that C4 may be chosen independently of r n , 
by Brownian scaling. 

Consider a reflecting Brownian motion X t = x + Wt + N t onD, starting from a point 
x G K\ U K 2 . Suppose that x + W t makes a double loop in an annulus B(x, r) \ B{x, 3r/4) 
for some r G (c 4 r n , r n /8), and then it makes a crossing from S(x + zr/2, r/16) to the ball 
B(x + i2r, r/16) within the convex hull of the two balls, before leaving B(x, r n /3), during 
a time interval [ti,^]- Suppose moreover, that L t2 — L tl < r/16. We will show that the 
two assumptions taken together yield a contradiction. The second assumption implies that 
\N t — N tl I < r/16 for all t G [ti, £2]- This implies that X will make more than one loop in 
B(x, 17r/16) \ B{x, llr/16) and then it will make a crossing from B(x + zr/2, r/8) to the 
ball B(x + i2r, r/8) within the convex hull of the two balls, before leaving B(x, r n /2). This 
is impossible because then X would make a closed loop around x within B(x,r n /2), and 
hence it would have to cross the boundary of D. We conclude that if the first assumption 
holds, then L t2 —L tl > r/16 > C4r n /16 = c$r n . Since the first event has probability greater 
than 1 — C3/2 and the process X starting from x E Ki U K 2 can hit 7 n _i before 7 n with 
probability greater than C3, the event that X hits 7 n _i before 7 n and L t2 — L tl > c^r n has 
probability greater than C3/2. This implies that reflecting Brownian motion in D starting 
from x G 7n-i/2 will hit 7 n _i before 7 n and Lj^ 1 > c$r n with probability greater than 
cq > 0. Hence, reflecting Brownian motion in D conditioned to hit 7 n _i before 7 n and 
starting from x G 7 n -i/2 will accumulate more than c^r n units of local time on dD n -i 
before hitting 7 n _i with probability greater than cq. 

Let A be the interior of D n _2 U -D n -i and let ip be a one-to-one conformal mapping 
of A onto a rectangle R = {a + ib : a± < a < a2,0< b < 1}, such that 7 n _2 is mapped 
onto the left side of R and 7 n is mapped onto the right side of R. Lemma 3.4 of [BCM] 
and a simple argument show that a 2 — ai is bounded above by a constant. Since the 
hyperbolic blocks are regular, there exists a point x G ip(D n _2) such that the harmonic 
measure of the upper part of R in ip(D n _2) is greater than c*, and the same is true for 
the lower part of the boundary. An analogous statement is true for ip(D n _i). All this 
easily implies that the distance of V>(7n-i) from the left and right sides of R is bounded 
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below by c-j = 07(0*) > 0. Let Ri = {a + ib : a\ + c 7 /2 < a < a 2 — C7/2, < b < 1}. 
By the Neumann boundary Harnack principle (Lemma 2.8), for any positive harmonic 
function h in R\ with Neumann boundary conditions on the upper and lower sides of R±, 
h(x) < cgh(z) for all tfj(-f n -i). This applies, in particular, to h y oi/j -1 . By conformal 

invariance, h y (x) < c%h y (z) for all i,zG 7n-i- 

Let <7(:r) be the harmonic function in D n _\ with Neumann boundary conditions on 
dD n _i fl <9-D, equal to 1 on 7„_i and equal to on j n . Reflecting Brownian motion in 
_D n _i conditioned to hit 7 n _i before 7 n is a g-transform of the unconditioned process. We 
have already proved that the g-process starting from x G 7 n -i/2 will accumulate more 
than c^r n units of local time on dD n -i before hitting 7 n _i with probability greater than 
cq. By the strong Markov property applied at the hitting time of 7 n _i/2, the same holds 
if the starting point belongs to 7 n . Without loss of generality, we may and do assume that 
hy(xo) = 1 for some xq G 7 n -i- Since < eg < g(x)/h y (x) < c\q < 00 for x G 7 n -i, 
an elementary argument shows that the /i y -process starting from x G 7 n will accumulate 
more than csr n units of local time on dD n -i before hitting 7 n _i with probability greater 
than cq. □ 



Proof of Theorem 3.2. (i) Let djy(x) = dist(x, Consider xq G -D. It is not hard to 
see that there exists zq G dD and a corresponding family of 7 n 's such that xq G D no for 
some no and dist(xo, dD n ) > cidr>(xo), where c\ G (0, 1) is a constant depending only on 
D. By the proof of Theorem 2.2 (see especially Lemmas 3.4 and 3.5) in [BCM], G{xq, ■ ) 
is bounded by 02k on for k < no — 1. Hence G(xo, • ) is bounded by 02^0 on D no -i. 
By the Harnack principle, it is bounded by c^uq on 8B(xq, cidr>(xo)/2), and since 



G(x ,x)=E x G(x ,X TB(xocidD(xo)/2) )j for x e D\B(x , Cl d D (x )/2), 
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the same bound holds on D\ B(xo, cido(xo)/2). We obtain, 

*t b 



E Xo L TBt = V" E Xo / l dDnndD (X t )dL t 

n=l - 70 

no-1 »T B „ °° /.Tb. 

= / . E Xo / lg DnndD (X t )dL t + N E Xo / l dDnndD (X t )dL t 

n=l - 70 n=n y ° 

ng — 1 „ oo „ 

= / C(xo, + / G(xo, x)cr((ia;) 

i JdD„ndD , „JdD„ndD 



n=1 JdD n ndD n=n<) JdD n ndD 

n — 1 



< |<9-D n n<9.D| sup G(x ,x)+ |<9-Dn n <9.D| sup G(x ,x) 
^ xedD n ndD xedD n ndD 

Tl() — 1 OO 

< \ dL> n n <9£>| c 2 n + |<9£>n n d£>| c 3 n 

n=l n=n 
oo 

< J^c 4 n|aD n naD|. 

n=l 

This is bounded by a constant independent of xq, by assumption (3.1). Hence we obtain 
sup xer) E x Lt b ^ < oo and, therefore, inf x(E £> u(x) = mf xe o E x exp(— £t b „ ) > 0. This 
means that the whole surface of D is active. 

(ii) Find a prime end € <9D and a family of 7 n 's such that (3.2) holds, that is, 

oo 
n=l 

Note that in the case of a simply connected domain D in R 2 , the Martin-Kuramochi 
boundary D* of D and the corresponding reflecting Brownian motion X* on D* can be 
realized as follows. Let cp be a conformal map from H = {a + bi : b > 0} to D and define 
D* to be the union of D and its prime ends. The map ip extends continuously to a one-to- 
one map from H to D*. Let Y be reflecting Brownian motion on H. Then (p(Y) is a time 
change of reflecting Brownian motion X* on D* . We will use this constructed reflecting 
Brownian motion X* in this proof. Recall that G{zq, x) is well defined for x E D* \ {zq} 
by the second paragraph of this section. For a > 0, define 

Va = {x E D* \ (B* U {^ }) : G(zq, x) = a} . 
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First, we claim that there exist positive integers mo, ni\ and a positive constant ao such 
that there is at least one D n , but at most m\ such sets, between r\ a and rj a+rno , for every 
a > ao. 

Recall that z* is the center of £?*. Let ip : H — > D be a one-to-one conformal mapping, 
such that </?(0) = zq and <p(i) = z*. Define h : D — > R by a(z) = — log |t/? _1 (z)|. Then a is 
harmonic in D with Neumann boundary conditions and a pole at zq. Let )j* = {i 6 D : 
h(x) = a}. Lemma 3.4 of [BCM] and conformal invariance easily imply that there exists 
an integer m such that for any a G R there is at least one D n between rf a and rf ~ . 

a-|-7TiQ 

It follows from the conformal invariance of the Green function that h\{z) = G{zq, <p(z)) is 
the Green function for reflecting Brownian motion in H starting from and killed upon 
hitting ip(B*). It is easy to see that h\{z) and — log \ z\ are comparable onHnjz : \z\ < r}, 
for some r > 0. This implies the existence of a positive integer mo and a constant ao > 
such that there is at least one D n between rj a and rj a+mo for every a > ao. From (3.1), 
the inequalities preceding (3.2) and (3.3) in [BCM] as well as Lemma 3.5 of [BCM], we 
see that there exists mi < oo such that there are at most mi sets D n between any rj a and 

Va+m ■ 

Let ctj be the sum of nr n restricted to integers n such that -D n -i lies between rj a and 
r] a+mo+ i, where a is of the form kmo +j. Every set D n _\ lies between r] a and rj a+rno+ i for 
some integer a, namely for the largest integer a such that i] a fl Q n -i = 0- This and (3.2) 
imply that J^™° a j = °°- We will assume without loss of generality that ao = oo. 

We define k(n) to be the integer k which maximizes kr^ among all fc's such that -Dfc-i 
lies between ?7( n _i) mo and r] nmo (we take the largest of the fc's with these properties if the 
above definition does not uniquely identify k(n)). If we restrict the sum in (3.2) to fc(n)'s, 
its value will be infinite, because there are at most mi sets D n _i between any r\ a and 
i] a+mo . By Lemma 3.5 of [BCM] and the comparability of — log \<p(z) \ and G(zo,z) for z 
in a neighborhood of zo, c\n < k(n) < c^n. 

By (3.1), the inequalities preceding (3.2) and (3.3) in [BCM] as well as Lemma 3.5 of 
[BCM], c\k < G(zo, x) < 02k for x E Dk for k > 1. Let j3 be the smallest integer multiple 
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of m greater than max{2, (c 2 /ci)}. We have 

oo f3 2j+m 

^2 Y nr fc(n) = °0 
j=l n =/3 2 J- 1 + m 

for m = or 1 and we will assume without loss of generality that we can take m = 0, i.e., 

oo 

Y nr k{n) = oo. (3.3) 

j=l n =p2j-i 

Let X* be reflecting Brownian motion on D* starting from some xq G D no , where no 
is large. Define 

Sj= mf{t>0:X* e Vf32j }, j>l, 
T(> n = inf{* > Sj : X* G rj nmo }, n>l, 
Ui> n ^ inf{t > T 3 k ,n : X; G 77 (n _ 1)mo }, n, k > 1, 
T^' n i f inf{* > : X* G Wo }, n,k>2, 

N£ = max{k : < ra > 1. 

In other words, iv~^ is the number of downcrossings of [(n — l)mo, nmo] by M t = G(-2o? ^t) 
between times Sj and Sj-i. This is of interest to us only for n such that [(n-l)mo, nmo] C 
[/? 2j— 2 , /3 2j ] . The process M is a continuous local martingale so it is a time-change of 
Brownian motion, until it hits 0. 

Consider a one-dimensional Brownian motion W starting from /3 2j and killed at the 
hitting time T of /3 2j_2 . It follows easily from the Ray-Knight Theorem that there is 
an event A with probability greater than pi > 0, such that on A, the local time Lj, 
accumulated by W at the level x before time T is greater than C4,9 2j for all x G (/3 2j_1 , /? 2j '). 
We will apply excursion theory to excursions of W from the set {nm : /3 2j_1 < nm < 
Z? 2 - 7 }. Given the local time {Lj,, x = nmo G [/3 2j_1 , /3 2j ]} and assuming the event A occurs, 
the distribution of the number of excursions going from nmo to (n — l)mo is minorized 
by a Poisson random variable with expectation K n > cs/^/mo = cq(3 2j . Conditional on 
{Lfi,£ = nmo G [/3 2j_1 , /3 2j ]}, these random variables are independent. Let T M (b) = 
in£{t >0:M t = b} and M J t = {M t+T M {p2j) ,t G [0, T M (/3 2 ^ 2 ) -T M {[5^)}}. Since M J t is a 
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time-change of W t , there exists an event A' with P Xo (A') > p 1 , such that on A', conditional 
on the local time of AF, the numbers of excursions of M- 7 between consecutive points of 
nmo < P 2 ''} are independent random variables minorized by independent 
Poisson random variables with means K n > cq(3 2j . 

Note that the processes AF are independent. We will now condition the process X* 
on the local times of AF's and the endpoints of excursions of X from {rj nmo , /3 2j_1 < 
nmo < P 2j }- 

Recall that k(n) is an integer such that -D n -i lies between 77( n _i) mo and rj nmo . An 
easy argument based on Lemma 3.3 shows that given endpoints of an excursion of X* 
going from r] nmo to 7/( n _i) mo , the amount of local time accumulated by the excursion on 
dD is greater than c^r^^ with probability greater than p 2 > 0. 

Let J n be the distribution of the local time accumulated by X* on the part of dD 
between ?7( n _i) mo and r] nmo , during the time interval (T M {f3 2j ),T M (P 2 ^ 2 )). We have 
shown that on an event Aj of probability greater than pi, J n is stochastically minorized 
by a random variable I n whose distribution is Poisson with mean greater than p2CQf3 2 i ■ 
c 7 r k(n)- Hence J n is minorized by a random variable I n which has mean A n greater than 
c 8 r fc(n)/3 2j and variance A n . Moreover, we can assume that the J n 's are independent given 
Aj. Hence, the local time accumulated by X between hitting of n^j and r\^-i, on the 
part of dD between these curves, is stochastically minorized by a random variable Hj 
such that on the event Aj, its mean is bounded below by J2j- f3 2 i- 1 <nm <i3 2: > c 8 r fe(n)/^ 2 ' 7 > 
Ylj- ^i- 1 <nm <)3 2 3 c 9 nr k(n) an d the variance is equal to its mean. It follows that Hj takes a 
value larger than bj = \ X^j /3 2 J- 1 <nmo</3 2 ^ c 9 nr fc(n) with probability greater than p 2 > 0. 
Since the M J, s are independent, we can assume that the Hj J s are independent. Let Aj be 
independent random variables with P(Aj = bj) = 1 — P(Aj = 0) = pi. Since the reflecting 
Brownian motion X* starting from xo G D nQ has to go through 7^ for j = no, no — 1, • • • , 1 

df 

before reaching 70 = dB*, the distribution of the local time accumulated by X* before 
hitting is minorized by the distribution of Y^=i^j- I n view of (3.3), Y^jLi^j = 00 ' 
and this easily implies that Yl'jLi^j = 00 > a - s - Hence, for any b G (0, 00), there is 
some n such that P (X]j=i A? > ^) > 1 — 1/^- This implies that for any x G D no , 
P XQ (LT Bsfc > 6) > 1 — 1/6. Therefore, inf x€ D E x exp(— Lt Ss( ) = and we see that part of 
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the surface of D is nearly inactive. 



□ 



Example 3.4. Our first example is very simple. Suppose that for some a > 1, 

D = {x = (xi, X2) : \x2\ < Xi and < X\ < 1} . 

The interesting range of the parameter is a > 1. We will show that if a G (1, 2) then the 

whole surface of D is active and when a > 2 then it is not. 

It is easy to see that it is sufficient to analyze only one boundary point, namely, (0, 0). 

We generate a corresponding system of hyperbolic blocks by letting 7 n 's be vertical cuts 

of the domain at distance 2~ k + j2~ ka from 0, for all j > such that 2~ k + j2~ ka < 
2 _ fc+ i _ 2 -fc^ for aU k > 2 

The number of hyperbolic blocks whose distance from is between 2~ k and 2~ k+1 is 
of order 2~ k ^ l ~ a \ Hence the blocks in this family have indices n of order Ylj<k 2~^ 1 ~ a ' ) ps 
2-fc(i-«) p er j me t er f eac h of these blocks is of order 2 _fca , so the contribution from 
these blocks to the sum in (3.1) is of order 2~ k ^-^ ■ 2- fc ( 1 " Q • 2~ ka = 2" fc ( 2 " a ). If a < 2 
then ^2 k>1 2 _fc ( 2_a ) < 00, so part (i) of Theorem 3.2 implies that the whole surface of D 
is active. 

The distance between 7 n and 7 n +i is comparable to the perimeter of _D n , so the same 
calculation as above shows that the sum in (3.2) is comparable to Ylk>i 2~ k ^ 2 ~ a ^ and this 
is infinite for a > 2. Therefore part of dD is nearly inactive when a > 2. 

The multidimensional version of this example will be discussed in Example 4.13. 

It is interesting to compare the above result with the semimartingale property of 
reflecting Brownian motion X in D starting from the tip = (0, 0). It is shown in DeB lassie 
and Toby [DT] that X starting from is a semimartingale if and only if a < 2. See also 
Theorem 3.1 (i) of Burdzy and Toby [BT] for a similar result. Fukushima and Tomisaki 
[FT] proved for domains in the shape of multidimensional cusps that reflecting Brownian 
motion starting from the cusp point is a semimartingale if a < 2. We will show in Remark 
4.14 below that it is not a semimartingale when a > 2. □ 
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Remark 3.5. In the definition of T>, it is required that \di,D\ be finite. One can of 
course relax this condition using localization. However if \di,D\ = oo (under whatever 
generalization one uses) and if u is a weak solution to (1.3)-(1.5) in the sense of (2.3) with 
< u < c, then m£ xe d L D u(x) = 0. For suppose otherwise, that is, there exists Co > 
such that u(x) > cq for every x £ d^D. Let g be a smooth function with compact support 
in IR d such that g = on and g = 1 on dD. Then by (2.3) we have 



This is impossible since the left hand side should be finite by the Cauchy-Schwarz inequal- 



Example 3.6. We will analyze a fractal domain which contains channels that become 
thinner at the same rate as the single channel in Example 3.4. In the present example, 
the distance between 7 n and 7 n +i is much smaller than the perimeter of D n for some n. 
Nevertheless, there is no gap between conditions (3.1) and (3.2) for this family of domains. 

Suppose that a > 0, (3 > 1 and let ak = Sj=i 2 _ ( J_1 ) a . Let S n be the family of all 
binary (zero-one) sequences of length n. We will write s = (si, S2, . . . , s n ) for s £ S n . For 
integer k > 1 and s £ <Sfe, we set b s = Y2j=i Let A* = [0, l] 2 , for k > 1 and s £ St 



A s = {(xi, x 2 ) £ M 2 : a k < x x < a k+1 , b s <x 2 <b s + 2 k/3 }, 

and let D be the connected component of the interior of A* U Ufc>i U s e<s fc ^ s ^ na ^ contains 
the open square (0, l) 2 (see Figure 3.2). (Note that when (3 > 2, the interior of A* U 
U*;>i Use5 fe is disconnected.) 
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The interesting range of parameters is/3>a>0. Ifa<l then \dD\ — oo, so part of 
the surface of D is nearly inactive by Remark 3.5. We will show that if 1 < a < (3 < 2a, 
then the whole surface of D is active and when (3 > 2a > 2, then part of the surface is 
nearly inactive. 

As in the case of Example 3.4, we will analyze only the family of hyperbolic blocks 
corresponding to a boundary point at the end of a channel. The analysis of other boundary 
points is straightforward but tedious so it is omitted. Fix a boundary point zq at the end 
of a channel, i.e., a point whose first coordinate is YlJLi 2~^~ 1 ' a . Let A k be the family of 
vertical lines K k , n = {{%, y) : x = cik + n2~ kf3 }, with n > 1 such that a k + n2~ kf3 < at+i- 
Let Ck be the family of these line segments in K kjn (~) D for K kjU £ Ak that separate zq 
from A*. Let {7^} be the relabelled family \J k C k . 

Let {s fc e S k ,k > 1} be the sequence such that the channel formed by the A Sfe 's 
approaches zq. The number of hyperbolic blocks D n defined by the 7 n 's needed to reach A Sk 
is of order X]j<fc 2 _JQ /2 _J/3 ps 2 k ^~ a K Consider a hyperbolic block D n which intersects 
A Sk The set D n may be either a square or it may contain a "tree" of thin channels. 
Consider first -D n 's that are squares. There are about 2 k ^~ a ^ such hyperbolic blocks, so 
they correspond to n in (3.1) of order ^7<fc ^ ' which is within a constant multiple 
of 2 k (P~ a \ The perimeter of any such D n is of order 2~ k @, so the total contribution of 
such IVs to (3.1) is of order 2 k{p -^2 k ^-^2- kp ^ 2 k( - f3 ~ 2a \ The series •£ k 2 k( -f 3 ~ 2 ^ is 
summable if and only if (3 < 2a. 

Next consider a D n which intersects A s with s £ S k and contains a side "tree" of thin 
channels. The length of its boundary is of order Ylj>k 2^~ fc ^2~ JQ pz 2~ ka . It corresponds 
to n in (3.1) of order 2 k ^~ a \ There are at most two such -D n 's for each A s , so their 
contribution to (3.1) is of order 2 k ^~ a ^2~ ka pz 2 fc (/ 3_2a ). Hence, the contribution of -D n 's 
with side channels is of the same order as the contribution of D n that have the square 
shape. We conclude that (3.1) holds if (3 < 2a. 

If (3 > 2a, then the contribution of the square -D n 's is enough to make the left hand 
side of (3.2) infinite, due to the estimates presented above. □ 

4. Higher dimensional domains. 
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This section is devoted to a family of multidimensional domains. The family may seem 
small, but it contains many examples that arise naturally in the context of the present 
paper and that of [BCM]. Before presenting the main result of this section, Theorem 4.3, 
we will state a definition and a technical assumption. 

Recall the family V from Definition 2.1. 

Definition 4.1. We will say that a domain D C M d , d > 3, belongs to the family V\ if 
\D\ < oo, D G V, and it satisfies the properties listed below. Recall the meaning of the 

df 

Lipschitz constant A from Definition 2.1. Let 70 = dB*. For every boundary point z G dD 
there exists a family of disjoint smooth (d — 1) -dimensional surfaces {'j n }n>i, such that 
7 n C -D, and the set D\-f n consists of two open connected components, fi n and Q! n . For 
every n, we have z G Q n , -B* C Q' n , and O n+ i C fi n . Let r n be the distance between 
7 n and 7 n +i- There exist ko < 00 and < cti, «2, . . . , 0.1 < 00, depending only on D, 
satisfying the following conditions. 

(i) For n > 0, a.\ < r n /r n+ i < and 7 n can be covered by at most ko balls of radius 

(ii) For every n > 0, there exists a curve r C D of length less than a^r n , connecting 7 n 
and 7n+i, whose distance from dD is greater than a^r n . 

(iii) For every n > and x G 7 n fl <9.D, there exists an orthonormal coordinate system 
CS with the property that dD n «6^n) is the graph of a Lipschitz function with 
constant A in CS. 

(iv) For every n > and x G 7 n fl there exists an orthonormal coordinate system 
CS with the property that dVt n fl B(x, aQr n ) is the graph of a Lipschitz function with 
Lipschitz constant a-j in CS, and the analogous statement is true for Q' n in place of 



We will write D n = Q n \ fl n +\. 

Note that it follows from part (ii) of Definition 4.1 that there is a constant ag > 
depending only on D G T>\ such that |£) n | > ctgr^ for every n > 0. 

Our proof of the second part of our main result in this section, Theorem 4.3, requires 
the following technical assumption, Condition 4.2. We will discuss ways of verifying this 
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assumption after the proof of Theorem 4.3. 



Condition 4.2. There exist < m < m\ < oo such that for any z G <9.D and 7 n 's as in 
Definition 4.1, if n > mi and xq G fl n , then sup x£7ii _ m[ G(xq, x) < inf x£7ii G(xq, x). Here 
G(x,y) is the Green function of reflecting Brownian motion in D killed upon hitting B*. 

Recall that we say that the whole surface of D is active if (1.6) holds. 

Theorem 4.3. Suppose that D C R d , d>3, is such that DeV 1 . 

(i) If for each boundary point z G dD, there exists a system of surfaces {^ n } as in 
Definition 4.1 such that 

n 

SUp \ dD ™ n 9D \ ^ < °°' ( 4J ) 

zedD n>l k=l 

then the whole surface of D is active. Here \dD n n dD\ denotes the (d — 1) -dimensional 
surface measure of dD n fl dD. 

(ii) Suppose now that Condition 4.2 holds. If there exists a boundary point z G dD 
and a family of surfaces {^ n } as in Definition 4.1, such that 

n 

E^E^-^oo, (4.2) 

n>l k=l 

then part of the surface of D is nearly inactive. 

The proof of the above theorem will be preceded by a few lemmas. Recall from Section 
1 that £?* C D is a fixed reference ball. In our proofs, Cj, kj, rrij and pj, j = 0, 1, . . ., will 
denote strictly positive and finite constants depending only on D. 

Lemma 4.4. Let D G T>\, zq G dD, and let {7™} and {r n } be defined relative to z$ as in 
Definition 4.1. There exist Ci, C2 G (0, 00), depending only on D, such that 

n n 

ciJ2rl- d <G(x,y)<c 2 J2r 2 k - d , 
k=i k=i 

for all n > 1, x G 7 n and y G 7 n +i- 
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Proof. Let {7^, n > 0} and {fi n , n > 0} be as in Definition 4.1, and recall that y) is 
the Green function for reflecting Brownian motion X* on D* killed upon hitting B*. As 
we observed in Section 2, X* = X on D U <9l-D. For /c > 0, let 



Gn fc (*,y) = G(x,y)-E x G(X T ^y) 



and note that Gn k (x,y) is the Green function for reflecting Brownian motion in Q k killed 
upon hitting 7/-. Note that since 70 = <9-B*, Ga (x,y) = G(x,y). 

It follows easily from Definition 4.1 that we can find points z k G D k , k > 0, and finite 
positive constants Co < ci depending only on D such that 

dist (z k , dD k ) > c r k and max {dist (z k , j k ), dist (z k , 7^+1)} < c x r k 

for every > 0. 

Let B k = B(z k ,cor k /2). Starting at any point in dB(z kl cor k /4), the expected time that 
Brownian motion spends in B k before hitting dD k is larger than c^r\. By the sup- 
port theorem for standard <i-dimensional Brownian motion, starting from any point in 
dB(z k , 3corfc/4), there is probability at least p± > (not depending on k) that the Brow- 
nian motion will hit the ball B(z k , cor k /4) before hitting dD k . So starting at such a point 
the expected time spent in B k before hitting dD k is at least "p\Civ\. This implies that 

/ Gn k (x,y)dx > pic 2 rl for every y e dB(z k , 3c r fc /4). 

Using the Harnack inequality and the fact that \B k \ = csrf, it follows that 

Gn k (z k , y) > c 4 r 2 k ~ d for every y e dB(z k , 3c r fc /4). 

By the Harnack and the Neumann boundary Harnack principle (Lemma 2.8), we have 

Gn k (x, y) > c 5 r 2 k ~ d , x G 7^+1 and y G 7 fc+2 . (4.3) 

On the other hand, starting in B k the expected amount of time reflecting Brownian 
motion X in D spends in B k before exiting the ball B(z k , 3cor k /4) is bounded by c§r\. By 
the support theorem for standard Brownian motion, there exists pi > such that starting 
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at any point in dB(z k , 3cor/j/4), there is probability at least P2 of hitting 7^ before hitting 
B k . So the number of crossing from B(zk, 3cor k /4) to B k by reflecting Brownian motion 
in fi fc killed upon hitting 7^ is majorized by a geometric random variable with mean 
l/p2- This implies that the expected amount of time spent in B k by X^ starting at any 
point in "fk+2 is at most cir\; that is 



/ C fife (x, y)dx < c 7 rl for y G 7 fc+2 . 

JB k 



Since \B k \ = c^rf, the Harnack inequality implies that 



G Qk {z k ,y) < c 8 r k 



2-d 



for every y G <y k +2- 



Again by the Harnack and the Neumann boundary Harnack inequality, we have 



Gn k (x, y) < c 9 rl d for every x G 7^+1 and y G 7^+2- (4.4) 

For k > 0, it follows from the strong Markov property that 

Gn k (x,y) = Gn k+1 {x,y) + E X G Qk (X T ^ +i ,y) for 

Consequently, for every a; G O^+i and y G Clk+2- 

G fifc (x, y) = Gn fe+1 (x, y) + E x [e w [G fife (X^ , y T v +2 )] ] , (4.5) 

where 7 is a reflecting Brownian motion in D independent of X and Z^ +2 is the first 

and note that 



hitting time of 7^+2 by Y. Let a k (x,y) = E x 
by (4.3) and (4.4), for x G flk+i and y G Qk+2, 



E, 



Gn fc (X T7fc+1 ,y T v +2 ) 



c 5 r z k < a k (x,y) < c 9 r 2 k d . 
Fix n > 1. By (4.5) and (4.6), for every x G 7 n and y G 7 n +i and < < n — 2, 

Gn k (x,y) - G Qk+1 (x,y) = a k (x,y). 
Adding these equations for < k < n — 2, we obtain 



n-2 



Gn (x,y) -Gn„_!(x,y) = ^a fc (x,y), 



fc=0 



(4.6) 
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or 

n-2 

G(x,y) = Gn (x,y) = Ga n ^{x,y) + ^a k {x,y). (4.7) 

k=0 

By (4.3) and (4.4), c 5 r^ < Gn^fay) < c 9 r^. This, (4.6) and (4.7) imply that 

n — 1 n—1 

c 5 J2rt- d <G(x,y)<c Q J2rl- d , 

k=0 fc=0 

for x G 7 n and y G 7 n +i- By Definition 4.1, «i < r n _i/r n < «2 and «i < tq/t-i < so 

n n 
k=l k=l 

for x G 7 n and y G 7 n +i- □ 

Lemma 4.5. For n > 3 and y G 7 n -3, iet x i— > /i y (a;) be the Poisson kernel with pole 
at y for reflecting Brownian motion in n _3 killed upon hitting 7 n _3. Let P| denote the 
distribution of the h y -transform of reflecting Brownian motion in D* killed upon hitting 
7 n _3, starting from x G n _ 3 . There exist Ci,pi > 0, depending only on D, such that 
P^(Lt 7ii _ 3 > cir n ) > pi for any collection of 7fc's as in Definition 4.1, any n > 3, x E "f n 
and y G 7 n _ 3 . 

Proof. All constants Cj that appear in this proof depend only on D. The conditions listed 
in Definition 4.1 imply existence of c 2 > and a point x G dD n _ 2 H <9L> such that (i) the 
distance from xq to 7 n -i U7„_2 is greater than 2c2r n , and (ii) there exists an orthonormal 
coordinate system CS Xo such that B(xq, C2r n ) fl 3D is the graph of a Lipschitz function 
with the Lipschitz constant A. Recall that A is the constant in the definition of V and, 
hence, in the definition of V\. We will assume that xq = in CS Xo and the positive part 
of the d-th coordinate axis intersects B(x , c 2 r n ) fl D. 

Let h(x) = P x (T 7ii _ 3 < T 7n ) and let xi be the intersection point of dB(xo, C2r n /2) and 
the positive part of the <i-th coordinate axis in CS Xo . It is easy to show, using Definition 
4.1 and a "Harnack chain of balls" argument, that h{x\) > C3 > 0. By the boundary 
Harnack principle (Lemma 2.8), we have h(x) > C4 > for all x G B(xq, 3c2r n /4) fl D. 
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Recall that X* is reflecting Brownian motion on the Martin-Kuramochi compactifi- 
cation D* of D, and X is the quasi-continuous projection of X* into D. As we noted in 
Section 2, X = X* on DUd L D. Let X* start from a point x G B(x , c 2 r n /2) n<9L> C <9l£>- 
It follows from Theorem 2.2 that Xt = x + Wt + Nt, where Wt is a ^-dimensional Brow- 
nian motion starting from and N t = J * n(X s )dL s is the singular push on the boundary 
8lD. Assume without loss of generality that A > 1. Let x 2 = (0, 0, . . . , 0, — c 2 r n /10), 
B 1 = B(x 2 , c 2 r n /(100A)), B 2 = B(0, c 2 r n /(100A)), and let C x be the convex hull of B 1 UB 2 . 
Consider the event A that the Brownian motion W hits B\ before leaving C\ in less than 
c^r^ units of time. By the support theorem and Brownian scaling, the probability of A 
is greater than p 2 > 0. Let T* = Tq B( ^ xq 3c2rn / 4 ) A c^r^. We will argue that if A occurs, 
then |iV T J > c 6 r n . To see this, first suppose that Tq B ^ 3car ^ < c 5 r^. Since A holds, 
W stays in Ci, and it follows that |WrJ < c 2 r n /5. Since — x| > c 2 r n /4, we have 
I-WtJ > c 2 r n /20. If Tq B( ^ xo 3car > and A holds, let to < be a time when 

W t() G B\. Since x + Bi is at a distance greater than c 2 r n /(100A) from L> and X to G -D, 
we must have \NtJ > c 2 r n /(100A). We see that with probability p 2 or greater, X accu- 
mulates at least cer n units of local time before leaving B(xq, 3c 2 r n /4). Since h(x) > C4 for 
all x G -B(x , c 2 r n ) fl D, X starting from any point x G B(x , c 2 r n /2) fl has a chance 
Ps > (depending only on D) of accumulating at least c§r n units of local time and hitting 
7 n _ 2 before hitting 7 n , by the strong Markov property applied at the time of hitting of 
dB(x ,3c 2 r n /4). 

Suppose that x\ G 7 n -i lies at least crr n units away from dD. By the support theorem 
for Brownian motion, the chance that reflecting Brownian motion X starting from x\ will 
hit B(xo, c 2 r n /2) fl dD before hitting any other part of dD U 7 n U 7 n _ 2 is greater than 
P4 > 0, depending only on D. By the strong Markov property applied at the hitting time 
of B(xo, c 2 r n ) fl dD, reflecting Brownian motion starting from x\ has a chance greater 
than p 5 > of accumulating at least c§r n units of local time inside B(xo, 3c 2 r n /4) and 
hitting 7 n - 2 before hitting 7 n . Hence, the /i-transform of X starting from x\ has a chance 
greater than p 5 > of accumulating at least c§r n units of local time inside B(xo, 3c 2 r n /4) 
before its lifetime. The boundary Harnack principle shows that the same is true for any 
x G 7n-i, except that the bound for the probability has to be replaced with a new value 
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P6 > 0. 

Now consider the /i y -transform of X starting from a point of 7 n . It must hit 7 n _i and 
then 7 n _2 on its way to 7^-3. The strong Markov property applied at the hitting times 
of 7n-i and 7^-2 and the claims proved so far show that an /i-process starting from any 
point in 7^ has a chance > of accumulating at least c§r n units of local time inside 
B(xq, 3c2r n /4) before its lifetime. By the Neumann boundary Harnack principle (Lemma 
2.8), there are positive constants C7 < c§ such that 



h(x) „ /i,y(:r) „ 
C 7T7 f — T / T — C 8T7 f forx,ze 7„_ 2 U7 n _i. 



A routine argument based on this observation allows us to extend the claim to the h y - 
process. □ 



Proof of Theorem 4.3. The main idea of this proof is the same as that of the proof of 
Theorem 3.2 but some details are different. 

(i) Consider x G D. It is not hard to see that there exists z G dD and a corre- 
sponding family of 7 n 's, as in Definition 4.1, such that xq G D no+ i for some no > 1 and 
dist(xo, dD) > cir no +i. For y G £l' no _i, by the strong Markov property of X, 



G(x ,y) = E 



Xq 



G(Xr Jno ,y) 



E 



Xq 



E, 



GlXn 



l fn 



-<Y 

-~in Q -l 



where Y is a reflecting Brownian motion in D* killed upon hitting dB* starting from y 



and independent of X, and TY is the first hitting time of 7^ by Y. Hence by Lemma 



4.4, G(xo,y) is bounded by C2 Ylk=i r k~ d f° r V e ^n -i- By the Harnack principle, y 1— > 
G(xo,y) is bounded by C3 X]fc=i r fc _d on dB(xo, cir no +i/2), and the maximum principle 
implies that the same bound holds onD\ -B(x , cir no+1 /2). We obtain, with cr denoting 
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the surface measure on d^D = \J^ =1 {dD n n 3D), 



l dLD (X t )dL t 
G(xo, x)a(dx) 

> / G(xo, x)cr((ia:) + / G{xq, x)a(dx) 

JdD n ndD JdD n ndD 



<d L D 
no — 1 



n=l J dD n ndD n=nQ ~„^ n , 
no — 1 oo 

< V |aD n n <9D| sup G(x ,x) + V* I^A* n <9L>| sup G(x ,x) 
n=1 xedD n ndD n=no xedD n ndD 

no — 1 n oo no 

< ]T |<9£> n n<9D|c 4 X>fe" d + E |0A»n0£>|c 3 X>j>- d 

n=l fc=l n=n fc=l 

oo n 

<^ C5 |<9L> n n<9L>| J>*- d . 

n=l fc=l 

This is bounded by a constant independent of x , by assumption (4.1). Hence we obtain 
sup^g^, E x Lt b < oo and, therefore, mf x& D u(x) = mf xe o exp(— £t b „ ) > 0. This 
means that the whole surface of D is active. 

(ii) Consider a point zq G dD and a corresponding family of 7 n 's satisfying (4.2). Let 
{x n ,n > 1} be a sequence in D that converges to zq. There is a subsequence {x nj ,j > 1} 
that converges to some Zq in D* , the Martin-Kuramochi compactification of D. Note 
that x I— > G(zq,x) is well defined on D* \ {zq} by the second paragraph of Section 3. In 
particular, 

G(zq,x)= lim G(x n . ,x) for x G D* \ {zq}. 

j^oo 

For a > 0, define 

77 a = {x G IT \ (B* U {z *}) : G(z *,x) = a}. 

Note that ?yo = <9-B*. Recall the definition of the integer tuq > 1 from Condition 4.2. Define 
for n > 1, 

a n = inf G(zq,£)- 

3? £E T3 71 777, Q 

By the Neumann boundary Harnack principle applied to the function x i— > G(zo,:r) on 
73nm ? there is a constant Co G (1, oo) depending only on D such that mf xej3nmo G(zq, x) > 

38 



coE x [G(z%,X T )] for every x G 73(n+i)m - F °r n > 0, 



a n+1 -a n = inf G^x)- inf G(zq,x) 

x £l3(n + l)m a;e73nm 



= c 



inf G 



S2 



3nm.Q 



inf 





C — 

Ifc73(™ + l)Tn 

2-d 



Go, (Xt , 



— C l r 3rem ' 



(4.8) 



where the last inequality is due to (4.4) and c\ > is a constant depending only on D. 
On the other hand, for x G 7 n , 



G(zo»aO = hm G(x nj ,x) = lim E x G(X T , 



So it follows from Lemma 4.4 that there are positive constants c 2 < c 3 depending only on 
D such that 

3iirao 3nmo 



.2-d 
k 



for every n > 1. 



(4.9) 



k=i u fc=i 

It follows from Definition 4.1 that for some C4 < 00, ?"n _c V r n-i < C4 for every n > 1. 
Hence, 



< 



E.inm 
k=l 



3nrao 2 — d 



r. 



\3nmo 



3(n-l)m 2-d — 



C 3 c 3 E 

fc=3(n— l)mo + l fc 



2-d 



C 2 



c 2 r 



2-d 

3(n — l)mo 



(4.10) 



<^(l + cf»)<oo. 

C2 



Let /3 = max{2, (3c 3 m /c 2 )(l + cf™ )} and n,- = inf{n : a n > ft}. By (4.10), 



ft < a ni < a n ,._i/3 < /?■ 



(4.11) 



Since by (4.2), 



3mQ — 1 1 oo 

ZEE E 

m=0 i=0 j=l {n:3nmo+m6(Ti4j_4 + 2i,n4j-2+2i]} 



3nmo+m 



r 



-d-l 
3nmo+m 



E 

fc=i 



2-d 

r fc = 00, 
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without loss of generality, we may and do assume the sum is infinite for m = and i = 1, 
i.e., 

oo 3nmo 

E E -ir™ E r *~ d = °°- ("- 12 > 

j=l {n:3nmo6(n4j-2,n4j]} &=1 

Let X* be reflecting Brownian motion on D* starting from some xq G D no , where no is 
large. Define 





= ini{t>0:X: e Van .}, 


3>1, 




= M{t> SjiX* erjaj, 


n>l, 


ur 


= mf{t > 4> n : X; G Van _ 


J, n,k>l 


k 


= inf{t > : X; G Van 


}, n, k > 2, 


N J 

1 'n 


= max{/c : U J k ,n < Sj-i}, 


n > 1. 



In words, iV^ is the number of downcrossings of [a n _i,a n ] by M t = G{z%,Xt) between 
times and Svj_i. This is of interest to us only for n such that [a n _i,a n ] C [a n -_ 1 , a n .]. 

The process M is a continuous local martingale so it is a time-change of Brownian 
motion, until it hits 0. 

Consider a one-dimensional Brownian motion Wt starting from a n4j and killed at 
the hitting time T of a n4j _ 4 . Note that by (4.11), a n4 . > /3 4j , a n4j ._ 4 < /3 4j ~ 3 , and 
a>n 4 j-2 > /3 4j ~ 2 . It follows from the Ray-Knight theorem that there is an event A with 
probability greater than pi > 0, such that on A, the local time accumulated by W at 
level x before time T is greater than c^fo for all x G (a n4 ._ 2 , a n4 .). We will apply excursion 
theory to excursions of W from the set {a n : < n < n4j}. Given the local time 

{Lf, : x = a n E {a n4j _ 2 , a n4j }} and assuming the event A occurs, the distribution of the 
number of excursions going from a n _i to a n is minor ized by a Poisson random variable 
with expectation Kj > C5 — . By (4.8), we have 

A' > c -iJ!L 

Conditional on {L^, x = a n G [a n4j .„ 2 , a n4j .]}, these random variables are independent. 

Let T M (b) = infO > : M t = b} and M/ = {M t+T M {an ^t G [0, T M (a n4j _ 4 ) - 
T M {a n4j )]}. Since M/ is a time-change of W t , there exists an event A' with P Xo (^4 / ) > pi, 
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such that on A', conditional on the local time of M J at {a n ,n4j_i < n < n^j}, the 
numbers of excursions of M J between consecutive points of {a n ,n 4j _i < n < n^j} are 
independent random variables minorized by independent Poisson random variables with 
means K n > ^ 2_ d ■ By (4.9) and (4.11), there is a constant cq > depending only on 
D such that 

V^ n 4j r 2-d 

Kj > c 6 — — d 

Note that the processes M J are independent. We will now consider the process X 
conditioned on the values of the following random elements: local times of the M J 's ac- 
cumulated at levels {a n ,n^j-i < n < n^j}, and the endpoints of excursions of X from 
{Va n ,n 4j -i <n< n 4j }. 

It follows from the definition of a n and Condition 4.2 that the following condition is 
satisfied 

There are at least three consecutive Dk 's between n an _ 1 and r\ an . (4.13) 

This and an easy argument based on Lemma 4.5 show that given endpoints of an 
excursion of X going from n an to n an _ 1 , the amount of local time accumulated by the 
excursion on dD is greater than c^rsnmo with probability greater than p 2 > 0. 

Let J n be the distribution of the local time accumulated by X on the part of dD 
between r] an _ 1 and r\ an during the time interval (T M (a n4j ),T M (a n4j _ 4 )). We have shown 
that on an event A, of probability greater than pi, J n is stochastically minorized by a 
random variable I n whose distribution is Poisson with mean greater than 

y*n 4 j r 2-d n 4j 

P2 c 6 fc 2-/ c 7 r 3nmo = p 2 c 6 c 7 ri~^ o ^ r 2 k ~ d '. 

r 3nm {) k=l 

Hence J n is minorized by a random variable with mean A n > cgr^^ Y^k=i r k~ d an< ^ 
variance A n . Moreover, we can assume that the 7 n 's are independent given A y Hence, the 
local time accumulated by X between the hitting of r\ an ^ and r\ an ^ 2 , on the part of dD 
between these surfaces, is stochastically minorized by a random variable Hj such that on 
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the event Aj, its mean is bounded below by 



E 



,2-d 
k 



{n:3nmo<E[ri4j -2 ,ri4j]} 



k=l 



and the variance equals its mean. It follows that Hj takes a value no less than 




E 



C 8 r 3nm E T < 



,2-d 
k 



{n:3nmoE[ri4j -2 ,ri4j]} 



with probability greater than ps > 0. 

Since the M J, s are independent, we can assume that the H^s are independent. Let 
Aj be independent random variables with P(Aj = bj) = 1 — P(Aj = 0) = ps- The 
distribution of the local time accumulated by reflecting Brownian motion starting from 
xq G Dn Ajo before hitting B* is minorized by the distribution of XljLi ^-j- m view of (4.12), 
Sj>i = °°5 an d this easily implies that Ylj>i ^-j = °°i a - s - Hence, for any b < oo, there 
is some jo such that P(X]j<j > b) > 1 — 1/b. This implies that P Xo (L Tb ^ > b) > 1 — 1/6 
for every x £ Dn 4jo ■ Therefore, inf xe £> E x [exp(— L Tb ^ )] = and we see that part of the 
surface of D is nearly inactive. □ 

We note that the only place where Condition 4.2 is used in the proof of Theorem 
4.3(h) is to prove (4.13). We will next discuss Condition 4.2 but first we need a lemma. 

Lemma 4.6. Let D e V\, z e dD and let {D n ,n > 0} and {7 n , n > 0} be as in Definition 
4.1. There exist a ,p > depending only on D such that the following holds, 
(i) For every n > 1 and every positive harmonic function h on the interior of D n _\ U D n , 
with Neumann boundary conditions on dD D -D n -i U D n , 



h(x) < a h(y) 



for every x, y e 7; 



(») 



For every n > 1, 



7n+l 



<T ln _ 1 ) > P0 



for every x G 7. 
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Proof, (i) Let h be a positive harmonic function on the interior of D n _i U _D n , with 
Neumann boundary conditions on dD fl D n -i U D n . It follows from parts (i) and (iv) of 
Definition 4.1 that there are k\ < oo and c\ > 0, depending only on D, such that 7™ can 
be covered by at most k± balls B(xk,cir n ). Moreover, for each one of these balls, either 
B{xk, 2cir n ) C D or B{xk, 2cir n ) HdD is the graph of a Lipschitz function. If x, y G 7 n and 
both points belong to one of balls B(xk 7 Cir n ), then there is a constant c 2 > depending 
only on D such that h(x) < c^hi^y), either by the usual Harnack principle (if the ball is 
inside D) or by the Neumann boundary Harnack principle proved in Lemma 2.8. It follows 
by a Harnack chain argument that h(x) < aoh(y) for any x, y G 7 n , where «o = c% • 

(ii) According to the definition of D G Pi, there is < 00 such that there ex- 
ists a "Harnack chain of balls" connecting 7 n and 7 n +i, that is, we can find a sequence 
S(xi,r),S(x 2 ,r-), . . .,B(x k ,r) in O n _i with k < k 2 , x x G 7 n , x fc G 7 n+ i and G 
B(xj-i,r/2) for j = 2, . . . , k. The existence of this "Harnack chain of balls" and the 
Harnack inequality easily imply that for some p± > depending only on D and some 
x G 7 n , 

Px(T 7n+1 <T 7n _J >pi. 

Applying part (i) of this lemma to harmonic function x 1— > P x (T 7ii+1 < T 7ti _ 1 ), we conclude 
that there is some po > depending only on D such that P x (T ln+1 < T 7ti _ 1 ) > po for 
every x G 7 n . □ 

Part (ii) of Theorem 4.3 has been proved under the assumption that Condition 4.2 
holds. Condition 4.2 seems to be difficult to verify in a direct way. We will state two 
other conditions, Conditions 4.7 and 4.8, that are easier to verify in examples. We will 
show that Condition 4.7 implies Condition 4.8 and Condition 4.8 implies Condition 4.2. In 
some examples, Condition 4.7 is the easiest condition to verify, but in some other examples 
Condition 4.8 holds even though Condition 4.7 does not. Lemma 4.11 below shows how 
one can verify Condition 4.7 in some examples. 

Condition 4.7. Let a and p be the constants in Lemma 4.6 and let T ln be the Grst 
hitting time of 7 n by reflecting Brownian motion in D. There exist < mo < mi < 00 
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such that for any z G dD and the 7 n 's as in Definition 4.1 corresponding to z, if n > n\\, 
p *( T 7n-m -i < T 7~+i) < a o 2 Po for every x G 7 n , 

and 

P x (T ln+1 < T 7ii _ mo _J < ao 2 p for every x G 7 n _ mo . 



Condition 4.8. There exist < m < mi < oo such that for any z G dD and the 
7 n 's as in Definition 4.1, the following is true for n > m\. Let A be the interior of 



U n - mo -i<k<n D k and p x (dy) = P x (T 7 „_ mo _ lU7 „ +1 G dy), for x G A. In other words, p x 
is harmonic measure on the set 7 n _ mo _i U 7 n +i inside A for Brownian motion reflected 
on dD. Then the Radon-Nikodym derivative d\i z jdp y < 1 on 7 n +i and d\x z jdp y > 1 on 

In— m — 1> f° r z ^ ln—m 

and y G 7 n . 

Lemma 4.9. Condition 4.7 implies Condition 4.8. 
Proof. Recall the constants «o and po from Lemma 4.6. 



Let A n be the interior of U n -m -i<fe<n D k and p x (dy) = Vx(X T „ /n _ m() _ lU „ /n+1 G dy), 
for x G A n . In other words, p x is harmonic measure in A n for Brownian motion reflected 
on dD. Fix a set C C 7 n +i- By Lemma 2.7, a; i— > p x (C) is a non-negative harmonic 
function of x G A n . By Lemma 4.6, 

Hx{C) < a p y (C) and p y (^ n+1 ) < a /M7n+i) 

for x, y G 7 n , and so 

Atx(C) < a* y A»x(7n+i) for x, y G 7n . (4.14) 

By Lemma 4.6 and Condition 4.7, for n > mi, / u x (7 n+ i) > po fo r ah x^Tn and ^2(7^+1) < 
«o 2 p for all z G 7n-m - Hence, 



Hz{ln+l) / Hx{ln+l) < «0 > 
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(4.15) 



for all z G 7 n _ mo , x G 7 n and n > mi. If reflecting Brownian motion in D starts from a 
point in 7 n _ mo , it has to hit 7 n before hitting 7 n +i. Hence, by the strong Markov property, 
(4.14), and (4.15), we obtain for n > mi, C C 7 n +i, z G r y n -m an d y G 7 r 



^(C) = / ^(C)P 2 (X TT _ mo _ iU7ri G cfa) 

^ 7n 



2 



= «o — 7 r/Wn+l) 

/V7n+l) 

< %(C)- 

Since C is an arbitrary subset of 7 n +i, the Radon-Nikodym derivative d\i z jd\i y < 1 on 
7 n+ i for z G 7 n - mo and y G 7 n with n > mi. Similarly, d[i z jd[i y > 1 on 7 n _ mo _i for 
z G 7n-m an d 2/ G 7 n with n > mi. Therefore Condition 4.8 is satisfied. □ 



Lemma 4.10. Condition 4.8 implies Condition 4.2. 

Proof. We will consider n > mi + 1. Suppose that xq G O n+ i and choose c\ so large that 
K = {x G D \ B* : G(xq,x) > Ci} C O n+ i. We will prove that sup x67ri mo G(xq,x) < 
inf x£7n G(x ,x). It will suffice to show that P X (T K < T B J < P y (T K < T B J for x G 
7n_ mo and y G 7 n , because G(;ro, x) = ciP x (T^ < T S J for x G D \ (B* U if). 

Our proof will use the technique of coupling. We will construct two reflecting Brownian 
motions in D on a common probability space, X starting from x G 7 n -m an d V starting 
from y G 7 n , such that {T% < Tjf } C {T£ < } almost surely, that is, 

P(r£<r£ and > t|"J = o. 

Let A n be the interior of Un-m -i<fc<n an d define 

^(cfy) = P x (Xr u edy) for x G A n . 

Given a; G 7 n -m an d y G 7 n , we will define some random variables on a common prob- 
ability space. Let rj X;n+ i be a random variable taking values in 7 n +i and having distri- 
bution fx x (dz)/ / u x (7 n +i). Let ?7y,x,n+i take values in 7 n +i with distribution (n y (dz) — 
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H x {dz)) / {n y {^ n+1 ) - // x (7„+i)), and let I x , n +i take values or 1, with P(I Xin+1 = 1) = 
/ u x (7 n+ i). Note that rj yjX}n+ i is well defined because, under Condition 4.8, n y (dz) > [i x {dz) 
on 7n+i- Similarly, let rj yjn - mo -i be a random variable taking values in 7 n _ mo _i and hav- 
ing distribution n y (dz) / ' ^, y {^ n - mo -i) on 7 n _ mo _i. Let r} x , v , n -m a -\ take values in 7 n _ mo _i 
and have distribution (/i x (dz)-ii y (dz))/(iJ, x (^ n - mo - 1 )-ii y (^ n - rno - 1 )). Let J y)n _ mo _i take 
values or 1, and assume that P(/ y)n _ mo _i = 1) = / u y (7 n _ mo _i). Due to Condition 4.8, 
we may and do assume that the J's are constructed so that I x ,n+i + I y ,n-m -i < 1? a.s., 
and we let I x>y = 1 — I x , n +i — I yjn - mo -i. Moreover, the ?y's are constructed so that they 
are independent, and independent of the J's. 

For x G A n and z G 7 n _ mo _i U 7 n +i, let denote the distribution of reflecting 
Brownian motion X in A n starting from x, conditioned on leaving A n \ {7 n _ mo _i U 7 n +i} 
through z. Let denote the distribution of a pair of processes (X,Y), such that the 
distribution of X is and the distribution of Y is Q^. The processes X and Y are 
defined on the same probability space but no further relationship such as independence is 
assumed. In particular, the two processes do not necessarily reach z at the same time. 

We will now define a distribution for a pair of processes (X, Y) starting from x, y G 
7 n _ m() U 7 n , such that either x = y or x G 7 n -m and y G 7 n . If x = y G 7 n -m„, we define 
= li for all t > 0, and the distribution of X is that of reflecting Brownian motion in 
.D, killed upon hitting 7 n . Similarly, if x = y G 7 n , define X t = Y t for all t > 0, and the 
distribution of X is that of reflecting Brownian motion in D, killed upon hitting 7 n _ mo , 
The most significant case is when x G 7 n -m an d 2/ G 7 n . In this case we let 

X Vx,n+lI X ,n+l ~(" Vy,n—mo — lIy,n—m(> — l ~"t~ Vx,y,n—m(, — l^x,y 

and 

^2/ — Vx,n+lI X ,n+l ~\~ f]y ,n— m () — 1 Iy,n— uiq — 1 ~f" T]y ,x ,n+ll x ,y • 

Note that by our construction, Z x and Z y have distributions ,u x and respectively. When 
Zy )n _ mo _i = 1 and Z x = Z y = ?? y)n _ mo _i = z G 7 n -m -i> we define the (conditional) 
distribution of (X,Y) to be Q x>y until the processes hit 7 n -m -i> an d then we "continue 
them as a single reflecting Brownian motion in D starting from z until it hits 7 n _ mo ." In 
other words, if X hits 7 n _ mo _i at time to and Y hits 7 n _ mo _i at time ti then ^t +t = ^tj+t 
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for t > 0. Similarly, when I x ,n+i = 1 an d Z x = Z y = i] Xjn+ i = z G 7 n +i, we define the 
(conditional) distribution of (X,Y) to be Q* until the processes hit 7 n +i, and then we 
continue them as a single reflecting Brownian motion in D starting from z until it hits 7 n . 
When J X)?y = 1 and Z x = z\ G 7 n _ mo _i and Z y = 22 G 7 n +i ; we let X have (conditional) 
distribution and then we continue it as reflecting Brownian motion in D starting from 
z\ until it hits 7 n _ mo , and we let Y be independent from X with conditional distribution 
Q^ 2 , and we continue it as reflecting Brownian motion in D starting from z^ until it hits 
7 n . We call the distribution of the processes constructed above P x , y - Note that under 
P x , y each one of the processes X and Y is a reflecting Brownian motion in D. Under P x ,y, 
the processes X and Y start from x, y G 7 n -m U 7n, i-e., Xo = x and Yo = y, they have 
random lifetimes ( x and ( Y , not necessarily equal, G 7 n -m U7„, and either 

= Yqy_ or G 7n-m an d ^c y - e 7n> The essential property of P x , y is that 

if X enters fl n +i before it is killed, then the part of the trajectory of X after the hitting 
time of O n+ i is a time shift of the trajectory of Y after its hitting time of Similarly, 
under P x , y , if Y enters fl' n _ mo _ 1 before it gets killed, then the part of the trajectory of 
Y after the hitting time of O n _ mo _i is a time shift of the trajectory of X after its hitting 
time of O n _ mo _i. 

We will use the distributions P x , y to construct processes X and Y which are de- 
fined on the whole time interval [0, 00). Suppose that x G 7 n -m an d y G 7 n and 
let (X 1 ^ 1 ) have distribution P x , y . Let (X 2 ,Y 2 ) have conditional distribution P X2 , V2 
given the event {X^ xl _ = X2 and Y^ Yl _ = 1/2}- We continue by induction. Given 
(X k ,Y k ), we let (X k+1 ,Y k+1 ) have conditional distribut ion P Xk+1: y k+1 given the event 
{X k xk _ = x k+1 and Y k Yk _ = y k +i}- It is easy to see that E& = 00 and Efc C^* = °°> 
a.s. Set C X ° = C Y ° = 0. For k > and t G [£o<,-<k C X \ E <;<fc+i C XJ ), define 

o<i<fc 

Similarly, for t G [ Eo<j<fc C Y ' , Eo<j<fc+i C^), define 

o<i<fc 
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It is straightforward to check that X and Y are reflecting Brownian motions in D and 
{T* < T£} C {Tj < TJJ. This proves that P X (T K < T B J < P y (T K < T B J for 
x G 7n-m an d 1/ 6 7n and, as we pointed out at the beginning of this proof, this implies 
that sup x€7ii _ mo G(x , x) < inf xe7 „ G(x , x). □ 

Recall A from Definition 2.1. 

Lemma 4.11. For any c\ and A there exists C2 such that the following holds. Suppose 
that for some n and we have \D^\ < cirf. for all n — — 1 < k < n. Then for all 



Px(T 7 „_ < T 7n+1 ) < C 2 : 



r 2-d 



Z^ii=n— m2 — 1 i 

and for all a; G 7 n -m 2 we have 

3-d 



n—m-2 — 1 



Px(T ln+1 < T ln _ m2 _ 1 ) < C2 n 2 _ (/ . 

2^ii=n— m,2 — 1 ^« 



Proof. We prove the second inequality, the first one being very similar. Write j for 
n — 1712 — 1. If a, b are integers with j < a < b < n, set L/^ = Ufc= a -^ fc ' define 

C a , b = inf{ / \Vf(x)\ 2 dx : / = on 7a and / = 1 on 7b+1 |, (4.16) 

and let R a ,b = C~\. C a ,b is called the conductance across U a ,b and R a ,b the resistance. 
Consider reflecting Brownian motion in U a ^ killed on hitting 7a and let G a ^(x,y) be the 
corresponding Green function. We use the fact that with respect to this process C a ^b is 
equal to the capacity of 76+ 1; see [FOT]. 

Using Definition 4.1 we can find a constant C3 independent of k and points Zk G Dk 
such that dist (zk, dDk) > c^rk- Let Bk be the ball of radius csrk/2 centered at Zk- 
Starting at any point that is a distance C3rfc/4 from Zk, the expected time that Brownian 
motion in spends in B& before hitting dD^ is larger than c^r\. By the support theorem 
for standard <i-dimensional Brownian motion, starting from any point that is a distance 
3c3rfc/4 from Zk, there is probability at least pi > (not depending on k) that the Brownian 
motion will hit the ball of radius C3rfc/4 about Zk before hitting dDk- So starting at such 
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a point the expected time spent in B k before hitting dD k is at least pic 4 rj?. Using the 
Harnack inequality and the fact that \B k \ = c^rf, it follows that G kjk (z k ,y) > CQr 2 ~ d if 
\y — Zk\ = 3c3rfc/4. By the Neumann boundary Harnack principle, 



G k ,k(z k ,y) > c 7 r 2 k a , ye lk+1 . (4.17) 

Consider reflecting Brownian motion in D k killed on hitting ^ k and let v k be the 
capacitary measure for 7fc+i- Then 



1 > Pz fe (T 7fe+1 < T lk ) = J G k , k (z k ,y)is k (dy) 



> 



crr 2 k d v k (lk+i) 
c 7 r 2 k ~ d C k , k . 



Therefore 



and 



i„rf— 2 



C kjk < c 7 r k 



R k ,k > c 7 r 2 k - d . (4.18) 

Next, if ai < (12 < (12 + 1 < as, let fx be the function on U ai ,a 2 a t which the infimum 
in (4.16) is attained and similarly f 2 the function on U a2+ x,a 3 - Let (3 = C a2+ i ;a3 /(C aija2 + 
C a2 +x,a 3 ) and define / on U ai ,a 3 by setting the restriction of / on U aijCl2 to be equal to 
(3 fx and the restriction of / on U a2 +x,a 3 to be equal to (3 + (1 — /5)/2- Then 

C ai , a3 < / \Vf(x)\ 2 dx = /3 2 f |VA| 2 + (1 - (3f I |V/ 2 | 2 



/3 C aua2 + (1 - I3Y C a2 + l,a 3 

C ai ,a 2 C a2 + x,a 3 



(4.19) 



Cai,a 2 + Ca 2 + l,a 3 

This is equivalent to 

Rai,a 3 ^ Rai,a 2 +-Ra 2 + l,a 3 - (4.20) 

By (4.18), (4.20), and induction, we obtain 



r 2-d 
' i i 



*=3 
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or 

Cj,n < =a 1 2 _ d . (4.21) 

Recall that Bj is the ball of radius c 3 rj/2 about z.,-. Starting in Bj the expected 
amount of time the process spends in Bj before exiting the ball of radius "Sc^rj/A about 
Zj is bounded by c%r 2 j. By the support theorem for standard Brownian motion, there 
exists P2 > such that starting at any point that is a distance 3csrj/4 from Zj, there is 
probability at least p 2 of hitting 7j before hitting dDj \"fj. A standard argument allows us 
to conclude that the expected amount of time spent in Bj starting at any point of Uj iTl is 
at most cgr'j. Since \Bj\ = cioTj , the Harnack inequality implies that Gj^ n (zj, y) < c\\r 2 ~ d 
if y G 7n+i- The Neumann boundary Harnack inequality then implies that 

Gj, n (x,y) < c 12 r 2 j~ d 7 x G 7j+i? V e 7n+i- ( 4 -22) 

Let z/ be the equilibrium measure for 7 n +i with respect to reflecting Brownian motion in 
Uj :n killed on hitting 7^. Combining (4.21) and (4.22), 

p x(T 7 „ +1 < T 7 .) = J G j>n (x, y)v{dy) < c 12 r)~ d C d ^ x G 7^+1. 

This proves the lemma. □ 



If the r n are comparable, then Lemma 4.11 implies Condition 4.7 for sufficiently large 

m . 

Remark 4.12. If D G X>i, then a "typical point" x G <9D has a neighborhood U C D 
such that <9-D n C/ is the graph of a Lipschitz function. The Green function satisfies 
G(x,y) < c\\x — y\ 2 ~ d , for x,y G U, where c\ depends only on the Lipschitz constant 
characterizing dD; to see this we flatten the boundary and reflect over a hyperplane as in 
the proof of Lemma 2.8, and then use the result of [LSW]. The upper estimate in Lemma 
4.4 follows from this immediately. 

Example 4.13. Our first example in this section is a multidimensional version of Example 
3.4. Suppose that d > 3 and for some a > 1, 

D = I a; = (xt, X2,...,x d ):Q<x 1 <l and x" > [x\ + . . . + a^) 1/2 j . 
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We will restrict the parameter range to a > 1. We will show that if a E (1, 2) then the 
whole surface of D is active and when a > 2 then part of the surface is nearly inactive. 

We will analyze only one boundary point, the origin, in view of Remark 4.12. We let 
the 7fc's be intersections of D with (d— l)-dimensional hyperplanes perpendicular to the first 
axis, at distances 2~ k + j2~ ka from 0, for all j > such that 2~ k + j2~ ka < 2~ k + 1 -2~ ka , 
for all k > 1. 

Note that for some c\ and any mo there exists m\ such that for any n > mi we have 
1/ci < Tj/rk < ci for all n < j, k < n + tuq. This and Lemma 4.11 easily imply that 
Condition 4.7 holds. 

The number of -D n 's whose distance from lies between 2~ k and 2~ k+1 is of order 
2 -fc(i-«). For D n 's in this range, Em=i r m" d ~ Ej<fc 2-J'( 1 - a )2-J Q ( 2 - d ) ss 2 - fc ( 1+Q! ( 1 - d )). 
The surface area, |9D n fl is of order 2~ fca ( d_1 ), so the contribution from these sets 
to the sum in (4.1) is of order 2~ k ^-^ ■ 2" fc ( 1 +«( 1 - d )) . 2~ ka ( d -^ = 2" fc ( 2 " a ). If a < 2 
then Efc>i 2 _fc ^ 2_a ) < oo, so part (i) of Theorem 4.3 implies that the whole surface of D 
is active. 

A similar calculation shows that the sum in (4.2) is comparable to Efc>i 2~ fc ( 2-a ) and 
this is infinite for a > 2. Hence, by Theorem 4.3 (ii), part of the surface of D is nearly 
inactive if a > 2. □ 



Remark 4.14. Fukushima and Tomisaki [FT] studied reflecting Brownian motion in 
unbounded cusps 

D=^x= (xt, x 2 , ■ ■ ■ ,x d ) : xt > and x% > {x\ + . . . + x^) 1/2 | 

and derived a Green function estimate (see Lemma 5.4 and 5.5 in [FT]). Their proof can be 
adapted to get the Green function upper bound estimate for reflecting Brownian motion 
in the truncated cusps D as defined in Example 4.13 and to show that for 1 < a < 2, 



/ G D \ B *(x,y)(r{dy) < oo. 

JdD 



>dD 

Thus this gives an alternative proof for the boundary of D to be active when 1 < a < 2. 
The main goal of the paper [FT] is to show that reflecting Brownian motion in D starting 
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from the cusp point = (0, • • • , 0) is a semimartingale when a < 2. Using Theorem 4.3(h) 
(and its proof) in this paper, we can settle the remaining case by showing that reflecting 
Brownian motion in D starting from is not a semimartingale when a > 2. Clearly, this 
is equivalent to the fact that reflecting Brownian motion in D starting from is not a 
semimartingale when a > 2. 

By Theorem 2.1 of [FT], reflecting Brownian motion X in D is a strong Feller process 
on D and thus can start from every point in D. Let a > 2. According to Example 4.13, 
part of dD is nearly inactive. By the proof of Theorem 4.3(h), 

lim P x {Ltb, > b) > 1 — - for every b > 0. 

Were X a semimartingale starting from 0, the Skorokhod decomposition for X 

X t = X + W t + [ n(X s )dL s fort>0 
Jo 

would hold under P x for every x G D. It follows from weak convergence and the second 
to the last display that 

P (L TSi = oo) = 1. 

This is a contradiction since P (Tb* < oo) > 0. Therefore X starting from cannot be a 
semimartingale. □ 



Example 4.15. This is a multidimensional analogue of Example 3.6. Suppose that a > 0, 
(3 > 1 and let = Y^=i 2~^~ l ^> a . Let S n be the family of all binary (zero-one) sequences 
of length n. We will write s = (si, s 2 , . . . , s n ) for s e <S n . For integer k > 1 and s e *Sfc, we 
set b s = J2j=i Sj 2 ~ j - Let A* = [0, l] d . For k > 1 and s G 5 fc let 

A s = {(xi, . . . , Xd) : a fe < xi < a fc+ i, ((x 2 - & s ) 2 + x\ + . . . + x^) 1/2 < / B (xi)}, 

where Ci2 _fc ^ < f s (xi) < 2~ kf3 and ci > does not depend on s. Assume that all functions 
f s are Lipschitz with the same Lipschitz constant. Let D be the connected component of 
the interior of A* U Ufc>i U s eS fe ^ s ^ n at contains the open box (0, l) d . 
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We restrict the range of parameters to (3 > a. Since we have assumed that a > and 
(3 > 1, the surface of D is finite and Remark 3.5 cannot be used to draw any conclusions. 

We will show that if a < (3 < 2a then the whole surface of D is active and when 
(3 > 2a then part of the surface is nearly inactive. 

We will analyze only a family of -D n 's corresponding to a boundary point at the end of a 
channel, in view of Remark 4.12. Fix a boundary point zq at the end of an infinite channel, 
i.e., a point whose first coordinate is YlTLi 2~^~^ a . Let Ak be the family of hyperplanes 
Kk,n = {(xii • • • 7 Xd) '■ x\ = ak + n2~ k/3 }, with n > 1 such that ak + n2~ kf3 < ak+i- Let 
Ck be the family of connected components of Kk, n H -D, for Kk, n £ Ak, which separate zq 
from A*. Let 7 n 's be the relabelled family \J k Ck- 

This assumption on the magnitude of / and Lemma 4.11 imply that Condition 4.7 
holds as long as the relevant -D n 's belong to the same A s . Lemmas 4.9 and 4.10 then prove 
that Condition 4.2, i.e., sup xejn _ m G{xq,x) < inf x€7ii G{xq, x), holds if 7 n _ mo and 7 n 
belong to the same A s . Then clearly Condition 4.2 holds in full generality if we replace 
mo with 2mo + 1. 

The number of -D n 's defined by the 7 n 's, needed to reach A s with s e Sk is of order 
J2j<k 2 ~ Ja / 2 ~ J(3 ~ 2 k ^~ a \ Consider a D n which intersects A s with s e S k . The set 
D n may either have diameter of order 2~ k " or it may contain a "tree" of thin channels. 
Consider first -D n 's that have diameters of order 2~ k/3 . There are 2 k ^~ a ^ such -D n 's, up 
to a constant, so for n in this range, £m=i r^ d w ^.< fc 2^/ 3 -«)2-J /3 ( 2 - d ) w 2 fc ^( d " 1 )- Q ). 
The surface area, fl is of order 2~ kf3( ^ d ~ 1 \ so the total contribution of such -D n 's 
to (4.1) is of order 2 k ^-^2 k ^ d -^-^2- k ^ d -^ w 2 fc ^" 2Q! ). The series E fc 2 fc ^" 2a ) is 
summable if and only if (3 < 2a. 

Next consider a _D n which intersects A s with s E Sk and contains a side "tree" 
of thin channels. Its surface area, \dD n fl dD\, is of order J2j>k 2 j ~ k 2~ ja 2~ j/3 ( d ~ 2 ^ fx 
2-fe(a+/3(d-2))_ There are at most two such -D n 's for each A s , so their contribution to (4.1) 
is of order 2 k ^ d -^-^2- ki * a+f3{d - 2 ^ fx 2 k ^~ 2a \ Hence, the contribution of D n 's with 
side channels is of the same order as the contribution of D n that have diameter of order 
2-*0. We conclude that (4.1) holds if (3 < 2a. 

If j3 > 2a then the contribution of D n : s with diameter of order 2~ k/3 is enough to 
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make the left hand side of (4.2) infinite, due to the estimates presented above. 



□ 



Example 4.16. We will analyze a multidimensional fractal domain vaguely resem- 
bling the von Koch snowflake, except that we will add barriers partly blocking the passage 
between the building blocks. Suppose that the dimension of the space is d > 3 and fix a 
parameter p e (0, 1/2). We will impose further restrictions on p below. For k > 0, let Ak 
be a finite family of open cubes with edge length with edges parallel to the axes, and 
satisfying the following properties. The family Aq consists of one cube A . The family A\ 
consists of 2d cubes which are disjoint from each other and are disjoint from A . One side 
of any cube in Ai lies on a side of A and these two sides of the two cubes have the same 
center. Now suppose that we have defined families Ak for k < n. Let A n be the union 
of all cubes in \Jk<n-^ k - Then A n +i is the maximal family of disjoint cubes that do not 
intersect A n and such that one side of each of these cubes lies on a side of a cube from the 
family A n , and has the same center. Let D* be the union of all cubes in |Jfc>o an< ^ no ^ e 
that this set is not connected because all cubes in this family are disjoint. We transform 
.D* into a connected open set by adding "passages" between cubes. Fix a parameter j3 > 1. 
For any pair of cubes which belong to A n -i and A n , and whose sides intersect and have 
a common center x, we add to D* the open ball B(x, 2p( n_1 ) /3 ). We let D be the union of 
D* and all such balls. Parts of the boundary of D are adjacent to D on both sides, and 
this is forbidden by Definition 4.1, strictly speaking. We could modify the domain D or 
even Definition 4.1 to cover this case, but that would be an unnecessary embellishment. 

We will determine for which values of p the surface area is finite because the example 
is not interesting if \dD\ = oo; in such a case a part of the surface is nearly inactive by 
Remark 3.5. The surface area of a cube with edge length p k is of order p k ^ d ~ r K The 
number of cubes in Ak is of order (2d — l) k . The total surface area of cubes in Ak is of 
order p k< ^ d ~ 1 \2d - l) k . The surface area of D is finite if ^fcP fc ^ _1 H 2 ^ _ l) fc < that 
is if p d ~ 1 (2d — 1) < 1. Hence, we are interested only in p less than (2d — l) _1 /( d_1 ). The 
function f(d) = (2d — l) _1 /( d_1 ) is increasing for d > 3 because, when we treat dasa real 
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argument, 

= 1 + (2d - 1) (log(2d - 1) - 1) 
(rf-l) 2 (2rf-l) d/(d - 1} 

for d > 3. We have /(3) = l/y/5, and lim^oo /(d) = 1. Hence, we can take p G 
(0, 1/2 A l/\/5) for any d. We will see that, as long as the surface area is finite, the value 
of p does not play any role in this example. 

As usual in our examples, we will analyze only a point z G dD that lies at the end of 
an "infinite" channel, i.e., such that any continuous path in D from the center z* of Aq to 
z must pass through at least one cube in every family At- Let T be a continuous path in 
D from z* to z that passes through a side of any cube in [J k>0 Ak at most once, and if it 
does so, then it passes through the center of that side. Let Zk be the intersection point of 
T and the side of the cube in Ak that is a part of a side of a cube in Ak-i- The curve V 
passes through all the z^s on its way from z* to z. 

For every Zk, let Ck be the family of all sets dB(zt, 2 J ) f)D, where j satisfies Ap^ k ~ 1 ^ < 
2i-i < 2 j+1 < p k ~ 1 /2. Note that each set dB(z k , 2 j ) n D contributes two sets to C k and 
each one of these sets is a spherical cap. Let C = \J k Ck and rename the elements of C as 
7 n , in the order in which they have to be passed on the way from z* to z within D. It is 
elementary to check that this family of 7 n 's satisfies the conditions listed in Definition 4.1. 

In this example, Condition 4.7 does not hold. We will argue that Condition 4.8 holds 
directly. Consider spheres dB(0, 2^'), dB(0, 2^+ fc »), dB(0, 2^+ ko + m °) and &B(0, 2^+ 2ko + mo ) 
and call them Si, Sz, S3 and S4. It is not very hard to prove that there exist large ko and 
m , such that P X (T S4 G A,T Si < T Sl ) < P y (T S4 G A, T Si < T Sl ), for A C S 4 , x G S 2 
and y G S3. We also have, for sufficiently large ko and mo, that P x (Ts 1 G A,Ts 1 < 
Ts 4 ) > P y (T Sl G A,T Sl < T S4 ), for A C Si, x G S 2 and y G S3, although the two 
claims are not symmetric and require somewhat different justification. By the reflection 
principle, for reflecting Brownian motion in D, P ' x {T ln+2kQ+mo G A,T ln+2kQ+mo < T 7 J > 
Py(T ln+2ko+mo G A,T Jn+2ko+mo < T ln ) for A C 7 n+2fco+mo , a; G 7 n+fco and y G 7 n+fco+mo , 
provided 7 n and 7n+2/so+m belong to the same family C k and lie on the same side of 
dD. We also have P X (T^ G A,T~ < T~ ) > P„(T T 6A,T 7 < . ) for 

A C 7 n , x G 7n+fc an d y e 7n+fc +m - If we fake only every /co-th element of the family 
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7™, this proves Condition 4.7 for 7 n 's which belong to the same family Ck and lie on the 
same side of dD. Hence, Lemma 4.10 proves Condition 4.2 for n restricted in such a way. 
However, this implies that Condition 4.2 holds for all n, with mo replaced by 2mo, for the 
same reason as in Example 4.13. 

The number of 7 n 's in C k is of order log(( / o fc ~ 1 /2)/(4p (A: ~ 1)/3 )) « k. If z\ G fi n+ i with 
sufficiently large n then by Lemma 4.4, the Green function G(z±, • ) can be bounded by 
ci J2j<k3P jl3< " 2 ~ d ^ — C2kp hf3 ( 2 ~ d ^ for x G D that lie between 7 n 's in C fc . The surface area 
of dD n fl dD corresponding to 7 n G Cfc is bounded by csp k( ^ d ~ 1 \ so the contribution of 
such ZVs to the sum in (4.1) is bounded by c 4 /c 2 / o A:/3(2 ~ d) /9 fc(d ~ 1) = c 4 /c 2 p fc(/3(2 ~ d)+d_1) . If 
/3 < (d - l)/(d - 2) then J2k fcV (/3(2_d)+d_1) < oo and Theorem 4.3 (i) implies that the 
whole surface of D is active. 

To find a lower bound for (4.2), we take into account only one D n corresponding to 
each family Ck, namely the one with the largest surface area. We obtain as a lower bound 
for (4.2) the quantity c^kp k ^ 2 ~^ p k ( d ~V = C5 fc p fc(/J(2-d)+d-i)_ \{ p > (d-l)/(d-2) then 
Sfc /cp fc< ^( 2 ~ d ) +d ) = oo so by Theorem 4.3 (ii), part of the surface of D is nearly inactive. 

5. Trap domains. 

The ideas developed in Section 4 allow us to prove a new result on "trap" domains 
introduced in [BCM]. The new result applies only to the class of domains V 1 presented in 
Definition 4.1 but that class contains some very natural examples of fractal domains, such 
as the multidimensional version of the von Koch snowflake presented in Example 4.16, 
that were not covered by theorems proved in [BCM] (see Example 5.2 below). There was 
a big gap between results on two-dimensional domains and higher dimensional domains 
in [BCM]. At first we thought the gap was purely technical in nature — complex analytic 
methods could not be used in higher dimensions. It turns out that the gap is in fact 
"real," in the sense that the multidimensional examples are considerably different from 
the two-dimensional examples — compare our Example 5.2 and Proposition 2.15 of [BCM]. 
Recall that £?* C D is a closed ball with positive radius and T B „ = inf{t > : X t G 
is the first hitting time of £?* by X. We say that D C d > 2, is a trap domain if 

sup E X T B = oo, (5.1) 

xeD 
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and otherwise D is called a non-trap domain. One can express (5.1) in a purely analytic 
way, namely, by saying that D is a trap domain if and only if 

sup / G(x, y)dy = oo 
xeD\B Jd\b 

See [BCM] for further discussion of basic properties of trap domains. 

Theorem 5.1. Consider a domain D E V\, D C IR d , > 3, with a finite volume. 

(i) If there exists a constant c < oo such that for each point z E dD, there is a system 
of surfaces {^ n , n > 0} as in Definition 4.1 satisfying 

oo n 

£lA,i£>r<<c, 

n=l fc=l 

then D is not a trap domain. 

(ii) If there exists a boundary point z E dD and a system of surfaces {^ n } as in 
Definition 4.1, such that 

n 

n k—1 

then D is a trap domain. 

Proof, (i) The proof is very similar to the proof of Theorem 4.3. Consider xq E D. It 
is not hard to see that there exists z E dD and a corresponding family of 7 n 's, as in 
Definition 4.1, such that xq E D no for some no and dist(xo, dD no ) > c\r n(i . By Lemma 
4.4, G(xq, •) is bounded by C2 Ylk=i r k~ d on ^n -i- By the Harnack principle, it is 
bounded by C3 Ylk=i r k~ d on 9B(xo,cir nQ /2), and the maximum principle implies that 
the same bound holds oaD\ B(xq, c\r no /2). For x E B(x G ,cir no /2) we have G(xq,x) < 
c 4 \x - x | 2 " d Efcli r fc~ d /( r no/ 2 ) 2 " d 5 by comparison with the Green function in R d . We 
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(5.2) 



(5.3) 



obtain, using the upper bound in Lemma 4.4 for n < no, 
/ G(xo,y)dy = y2 G(x ,y)dy 

= / G ( x o^y) d y+ / G(x ,y)dy 

l<n<n J °™ n>n J D n\B(x , Cl r nQ /2) 

+ / G(x ,y)dy 

J B(x n ,CTr„ n /2) 



B(x ,c 1 r n(j /2) 

n no no 

2-d 



l<n<no fc=l n>no fc=l fc=l 



n 



r 2-d 

n>l k=l 



This is bounded by a constant independent of xo, by assumption (5.3). The theorem 
follows in view of (5.2). 

(ii) A calculation similar to that in part (i), based on the lower bound in Lemma 4.9, 
easily implies part (ii) of the theorem. □ 



We would like to emphasize that part (ii) of Theorem 5.1 is much easier to prove 
than part (ii) of Theorem 4.3. This is because all we have to show is that the function 
x I— > E x [TgJ is unbounded. In the proof of Theorem 4.3 (ii) we had to prove that the 
random variable Lt b ^ for reflecting Brownian motion X* starting from x converges to 
infinity in distribution as x approaches a boundary point zq G dD. 

Example 5.2. Recall the domain D and notation from Example 4.15. Recall that 
the number of 7 n 's in Ck is of order log((p fc_1 /2)/(4p( fc_1 ) /3 )) « k. If z\ G O n+ i with 
sufficiently large n then by Lemma 4.4, the Green function G(zi, ■) can be bounded by 
c i J2j<k3P J0(2 ~ d) < c 2 /c/9 fc/3(2_d) for x G D that lie between 7n 's in C k . The volume of a 
D n corresponding to a 7„ G Ck is bounded by c^p kd , so the contribution of such D n : s to 
the sum in (5.3) is bounded by c 4 k 2 p k ^ 2 ~ d ^ p kd = c 4 fc 2 /9 fc(/3(2 " d)+d) - If /? < d/(d-2) then 
Y2k p fc (^( 2_d ) +d ) < oo and Theorem 5.1 (i) implies that D is not a trap domain. 

To find a lower bound for (5.4), we take into account only one D n corresponding to each 
family Ck, namely the one with the largest volume. We obtain as a lower bound for (5.4) the 
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quantity c 5 kp k ^ 2 - d ^ p kd = c 5 kp k ^ 2 - d ^ +d l If (3 > d/(d-2) then £ fe kp k ^ 2 - d ^ = oo 
so by Theorem 5.1 (ii), D is a trap domain. □ 



We will use the above example to compare Theorem 5.1 to a result about multidi- 
mensional trap domains proved in [BCM]. The result in [BCM] was based on the notion 
of a J a -domain, used by Maz'ja in his book on Sobolev spaces [Maz]. Here is an informal 
definition of a J a -domain (see [Maz] or [BCM] for the rigorous definition). We say that 
D is a J a domain if for every smooth (d — l)-dimensional surface A which divides D into 
two connected components D± and -D^, we have min (\D£\, \D£\) a < ci|A|, where a < oo 
depends only on D (here |A| is the (d — 1) -dimensional surface area). Theorem 2.4 of 
[BCM] implies that if D is a J a domain with a < 1 then D is not a trap domain, and there 
exists a trap domain D e J\. 

Roughly speaking, one can determine whether the domain D of Example 5.2 belongs 
to J Q with a given a by comparing the surface area of the opening between cubes in Ah 
and Ah-i to the volume of the cubes in \Jj >k Aj. The surface area is of order p fc £( d_1 ) 
and the volume is of order p kd , so D is a J a domain with a < 1 if (3 < d/(d — 1). Hence, 
for the family of domains in Example 5.2, Theorem 2.4 of [BCM] shows that D is not a 
trap domain if (3 < d/(d — 1), while Theorem 5.1 of this paper shows that this holds for 
all j3 < d/(d — 2), and in addition it shows that this result is sharp. The gap between 
the power of the two approaches is not as striking in dimensions d > 3 as it is in the 
2-dimensional case, discussed in Proposition 2.15 of [BCM]. 
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